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The observational status of inflation after the Planck 2013 and 2015 results and the 
BICEP2/Keck Array and Planck joint analysis is discussed. These pedagogical lecture notes 
are intended to serve as a technical guide filling the gap between the theoretical articles on 
inflation and the experimental works on astrophysical and cosmological data. After a short 
discussion of the central tenets at the basis of inflation (negative self-gravitating pressure) 
and its experimental verifications, it reviews how the most recent Cosmic Microwave 
Background (CMB) anisotropy measurements constrain cosmic inflation. The fact that 
vanilla inflationary models are, so far, preferred by the observations is discussed and the 
reason why plateau-like potential versions of inflation are favored within this subclass of 
scenarios is explained. Finally, how well the future measurements, in particular of i?-Mode 
CMB polarization or primordial gravity waves, will help to improve our knowledge about 
inflation is also investigated. 
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I. INTRODUCTION 

With the release of the Planck data 2013 [1-3] and 2015 [4—6], and the recent 
BICEP2/Keck Array and Planck joint analysis [7], the theory of cosmic inflation [8— 
14] has acquired a new status. Several of its predictions such as spatial flatness of our 
Universe, the presence of Doppler peaks in the Cosmic Microwave Background (CMB) 
multipole moments, almost scale invariant power spectrum for density perturbations 
have been definitively confirmed by the recent CMB anisotropy measurements. That 
makes inflation a predictive and verified theory of the early Universe. 

In fact, another remarkable outcome of the Planck data is that they also allow us to 
identify which version of inflation is most likely to have been realized in Nature [15—17]. 
As is well-known, inflation comes in different flavors but these different scenarios make 
different predictions and, thus, one can, at least in principle, distinguish among them. 
The fact that the primordial fluctuations are adiabatic and Gaussian to a relatively 
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high degree of accuracy [2, 3] is an important indication that we probably deal with 
single-field slow-roll inflation (with standard kinetic term), the simplest bnt non-trivial 
model of inflation. Of course, the final word has not yet been spoken since many non¬ 
vanilla inflationary scenarios are still compatible with the data. But, presently, they 
are just not needed in order to explain CMB measurements even if this situation could 
change in the future. 

The fact that we now have high accuracy CMB data at our disposal also allows us 
to detect the “fine structure” of inflation and to constrain the shape of the inflaton 
potential. Here again, the Planck data have provided precious information. We 
now know that the potential is of the plateau type and that simple monomials are 
disfavored [15-18]. Moreover, we now start probing the reheating epoch [19, 20]. 
Reheating is the epoch, after inflation and before the radiation dominated era of the 
standard hot Big bang phase, where the inflaton field decays and where all matter we 
see around us was produced [21-24]. It is therefore of major conceptual importance. 
And Planck 2013 and 2015 data put non trivial constraints on the physical processes 
that took place at that time [19, 20, 25—29]. 

The goal of these lectures, given at the second Jose Plinio Baptista school on 
Cosmology held in Pedra Azul (Brazil) in March 2014, is to review how the above 
conclusions can be established. Many reviews on inflation can be found in the 
literature [30-33] and there are technical papers reporting the astrophysical and 
cosmological observations, snch as the Planck papers [1-3]. But, in between, few things 
can be found and the present article aims at filling this gap. In some sense, it can be 
viewed as a technical guide which, from a reasonable prior knowledge of inflation, 
permits a detailed understanding of the implications for inflation of the recent high 
accnracy CMB data. 

These lecture notes are also written at a special time: the Planck 2013 and 2015 
data [1-6] have been released and their conseqnences (in fact, mainly the consequences 
of Planck 2013) already analyzed in several works. Moreover, very recently, a joint 
analysis made by the BICEP2/Keck Array team and the Planck collaboration [7] 
has been published showing that the BICEP2 detection of B-mode CMB polarization 
announced in Ref. [34] is mainly due to dust and cannot be attributed to primordial 
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gravity waves produced during inflation. At the time of writing, the Planck 2015 
scientific products (in particular, the likelihood) are expected to be delivered in June 
2015 only. This means that reproducing or extending the Planck 2015 analysis is not 
yet possible. However, from what is already known, the Planck 2015 results are in 
good agreement with Planck 2013. Therefore, the conclusions discussed in the present 
article (model comparison, constraints on reheating etc ...) will most likely remain 
valid for the second release of the Planck data. Whenever available, we quote the 
values obtained by Planck 2015. 

These lectures are also related to the lectures given by C. Byrnes on Non- 
Gaussianities [35] and by D. Wands on CMB physics. Hopefully, these three reviews 
should provide the reader with a rather complete overview of modern primordial 
Cosmology and its observational implications. In particular. Ref. [35] reviews how 
Non-Gaussianities are produced in non-vanilla inflationary models while, here, we 
restrict ourselves to simple scenarios for which Non-Gaussianities are very small. The 
two lectures are therefore complementary. The lecture notes by D. Wands explain in 
details how CMB anisotropies are generated while, here, we just take it as a known 
fact (see also the recent review [36]). Therefore, the present article and the one on 
CMB physics are also complementary. 

These lecture notes are organized as follows. In the next section. Sec. H, we present 
general considerations on inflation. Rather than discussing inflation in detail, which 
can be found in many review articles, we just give the basics and choose to focus 
on the fundamental principles at the basis of the inflationary mechanism and its 
experimental justifications. In Sec. HI, we discuss how inflation can be realized in 
practice. In particular, in Sec. HI A, we review how inflation can be embedded in 
high energy physics. Recently, alternative parametrizations have been considered and 
we discuss them in Sec. HIB. In Sec. HIC, we also review how the reheating phase 
can be described. Then, in Sec. IV, we discuss the theory of inflationary cosmological 
perturbations of quantum-mechanical origin. This part of the inflationary scenario 
is especially important because this is how one can relate theoretical predictions to 
astrophysical observations. In Sec. IV A, we present the calculation of the two-point 
correlation functions, or power spectra, for scalar and tensor perturbations in the 
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slow-roll approximation. In Sec. IVB, we review the calculation of the three-point 
correlation function, or bispectrum, and in Sec. IV C, the calculation of the four- 
point correlation function, or tri-spectrum. All these considerations are made in the 
slow-roll approximation and for single-field models with minimal kinetic terms. In 
Sec. IV D, we discuss the isocurvature perturbations and how they can be produced 
in the framework of inflation. In Sec. V, we use the tools introduced before and 
compare the inflation predictions to the high accuracy CMB Planck data. In Sec. V A, 
we consider the measurements of spatial curvature, in Sec. VB the measurements of 
isocurvature perturbations and, in Sec. V C, those of Non-Gaussianities. Since these 
data indicate that single field models are preferred, we then focus on this class of 
scenarios. In Sec. VD, we give the constraints on the slow-roll parameters and on the 
derived power-law parameters, such as the spectral index, the running or the tensor- 
to-scalar ratio. We also discuss the implications of the recent joint analysis made by 
the BICEP2/Keck Array team and the Planck collaboration. In Sec. VE, we carry out 
a Bayesian analysis to do model comparison and determine what are the best models 
of inflation. In Sec. VF, we present the constraints on reheating that can be inferred 
from the Planck data. Finally, in Sec. VI, we recap our main results and discuss which 
lesson can be drawn for our understanding of inflation and primordial cosmology. 


II. GENERAL CONSIDERATIONS ON INFLATION 


The motivations for introducing a phase of inflation, i.e. a phase of accelerated 
expansion, are well-known; postulating a > 0 (a is the Friedmann-Lemaitre-Roberston- 
Walker -FLRW- scale factor) allows us to avoid the puzzles of the standard hot Big 
Bang theory (for a detailed discussion of these issues, see Refs. [30, 37, 38]). If gravity is 
described by General Relativity (GR), then, in a homogeneous and isotropic Universe, 
the equations of motion are given by 
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where pi and pi are respectively the energy density and the pressure of the fluid 
“i”. In the standard model of Cosmology, we have indeed a collection of different 
fluids, pressure-less matter (made of baryons and cold dark matter), radiation (made 
of photons and neutrinos) and dark energy. These different types of matter source 
the Einstein equations and control the dynamics of the expansion. Notice that the 
expansion rate of the Universe is given by the Hubble parameter which, according to 
the above equations, is defined hy H = a/a where a dot means a derivative with respect 
to cosmic time. The quantity Mpi is the reduced Planck mass and, in the following, 
we will also use the quantity k = 1/M|j = being the Newton constant. 

Finally, the quantity fC, that can always be normalized to 0 or ±1, represents the 
curvature of the spatial sections. Notice that one can also define an effective curvature 
energy density by pcurv = — 3/C/(Ka^) such that the Friedmann equation takes the form 
= (k/3) Y^iPi + (K/3)pcurv Defining Dj = p/pcri and fl/c = Pcuw/Pcri, where the 
critical energy density is pcri = /k, the Friedmann equation can be rewritten as 

Yli ^i + ^K. = 1 - 


Let us now discuss under which physical conditions inflation can be obtained. The 
above equations can be combined and lead to the following formula which relates the 
acceleration of the expansion to the matter content of the Universe 


a 
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This immediately implies that, in order to have inflation, the pressure must be negative, 
i.e. p < —p/3 where p and p are defined in Eqs. (1.1) and (1.2). Having realized that 
we need a negative pressure, the next question is of course which kind of matter can 
possess this property and this will be the subject of the two next sections. Of course, 
as is well-known, we will see that scalar fields are ideal candidates. 

But before starting this discussion, it is interesting to notice that inflation is a 
genuine relativistic effect since it involves the term 3p in the above equation (1.3), 
which is absent in Newtonian physics. Indeed, let us consider a sphere of radius R{t) 
and of uniform density p. A galaxy of mass m, located at the edge of the sphere, 
feels a gravitational field G that can be simply evaluated by means of the Gauss’s law, 
f G ■ dS = AttG^M, where M is the mass of the sphere. This gives G = G^M/B?. As 
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a consequence, the acceleration of the galaxy can be written as 




(1.4) 


or 


R _ 4:7rG^ _ p 
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(1.5) 


where we have used M = AirpR?/^. This equation is similar to Eq. (1.3) except that 
the term 3p is not present. The physical reason behind the presence of this term is 
deeply rooted in the fundamental principles of GR: in GR, every form of energy weighs, 
including pressure. 

The term 3p is so important for inflation that it is interesting to ask whether 
it plays a role in other physical situations and if its appearance has been tested 
experimentally and/or observationally. This is a difficult question since, in ordinary 
cases, the contribution of pressure is usually negligible, p <C p. In fact, four situations 
where a gravitating pressure is important can be identified: inffation, dark energy but 
in some sense this is the same as inffation, neutron stars and Big Bang Nucleosynthesis 
(BBN). In particular, it is interesting to see what can be said about the 3p terms in 
the last two examples. 

Let us start with the internal structure of a neutron star [39]. As is well-known, 
it is controlled by the Tohnan-Oppenheimer-Volkoff equations that can be obtained in 
the following way. The metric for a static and spherically symmetric solution can be 
written as 


ds2 = + e^^dr^ + (d^^ + sin^ Odip^) , (1.6) 

where t is time, r a radial coordinate and 0 and p angular coordinates. The quantities 
<I> and A are functions of r only. Matter is assumed to be described by a perfect fluid, 
the stress energy tensor of which can be expressed as 

Tfj.u = ip + p)Ufj,Uu+pgfj.u, (1-7) 

where g^i, is the metric tensor and the normalized 4-velocity reads = (—e^, 0, 0, O). 
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Then, the time-time and r — r component of the Einstein equations read 
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On the other hand, energy conservation, = 0, more precisely its radial 

component, implies that 


dp 

dr 


ip+P) 


d$ 

dr 


( 1 . 10 ) 


The other components lead to the fact that p does not depend on time, 9 or p, that is 
to say p = p{r). If we now define the gravitational mass m(r) by 


G^rn{r) = ^ (l - e , 


( 1 . 11 ) 


then Eq. (1.9) implies that 

= 47rp(r)r^. (1-12) 

dr 

Introducing the expression of the mass (1-11) in Eq. (1.8) in order to express dd>/dr 
and, then, inserting the corresponding expression in the conservation equation (1.10) 
leads to 


dp 

dr 


-(P + P) 




r2(l — 2mG^/r) 


m{r) + 3p(r) 



(1.13) 


The important point in this formula is that the term 3p participates to this expression. 
This means that self-gravity of pressure affects the internal structure of the neutron 
stars. In practice, in order to calculate this internal structure, one has first to choose 
an equation of state p = p(p). Once this is done, one can integrate the two following 
equations 


dr 


dp dp 
dp dr ’ 


dm 

dr 


47rp(r)r^, 


(1.14) 


the last equation being nothing but Eq. (1.12). This leads to the functions p(r) and 
m(r). The radius of the star, Rstar, is defined by p (i?star) = 0 and its mass is given by 
Mstar = fn (i?star)- One Can then plot the mass-radius relation Mgtar(.Rstar)- Of course. 
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FIG. 1. Mass-radius relations of neutron stars for different equations of state (“standard” in 
the left panel, more “exotic” in the right panel). Black curves correspond to the standard GR 
calculation while red curves represent the case where self-gravity of pressure is absent. Figure 
taken from Ref. [39]. 


one obtains different mass-radius relations for different equations of state. Let us also 
notice that, at fixed equation of state, one obtains a curve, and not a unique prediction, 
because one needs to specify p{r = 0) to be able to integrate the above equations. One 
thus has a family of points parametrized by p (r = 0). Several examples are displayed 
in Fig. 1 (black lines). 

The fact that the structure of a neutron star depends on the general relativistic 
3p term opens the possibility to experimentally test it. In order to do so, the idea 
of Ref. [39] is to study an ad-hoc modification of the Tolman-Openheimer-Volkoff 
equation such that 


^ = -( ^ 

dr P P 7-2(1 — 2mGf^/r) 


m{r) + 3xp{r) ( -vrr^ 
o 


(1.15) 


where x is a new, phenomenological, parameter introduced by hand. The term 3p 
weighs normally when X = 1 and does not weigh at all when x = 0- Notice that x = 0 
is not the Newtonian limit because there are other relativistic terms in Eq. (1.15) (for 
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instance 1 — 2mG^jr at the denominator). So the idea is now to re-derive the mass- 
radius relation for neutron stars and to see the influence of a parameter x 7^ 1; the 
hope being to be able to put constraints on x from astronomical observations. The 
results are shown in Fig. 1. The fact that red curves (namely those obtained with 
X = 0) are different from the black ones (those obtained in the standard GR case) 
confirms that the 3p term has a significant influence of the mass-radius relation. 

However, as shown in Fig. 2, the fact that the equation of state is not known 
accurately completely blurs the effect. Indeed, one sees that the corresponding 
uncertainty is typically of the same order of the effect we try to detect. Therefore, 
the conclusion is that, although it is true that self-gravity is crucial in order to predict 
correctly their mass-radius relation, at least for the moment, neutron stars cannot be 
used to experimentally test the 3p term. 

Let us now turn to the other possibility, namely BBN [40]. Since BBN takes place 
during the radiation dominated era for which p = p/3, it is clear that the 3p term 
should have an important impact on BBN. In order test the influence of the 3p term, 
we follow the same strategy as for neutron stars and introduce an ad-hoc modification 
of GR characterized by the y parameter, namely 

This equation should be compared to Eq. (1.3). In order to derive the Friedmann 
equation, we need another equation and we can use the first law of thermodynamics 
for an adiabatic expansion, namely d{a^p) = —pd(a^), written for a co-moving volume 
or, equivalently, p + 3a{p + p)/a = 0. Then, noticing that d/a = l/(2da)d(d^)/dt and 
using the conservation equation, it is straightforward to derive the following relation 

d (a^) = [(1 - 3x) pada - ya^dp] . (1.17) 

If X = 1, it is easy to check that 
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FIG. 2. Mass radius relations for different equations of state and associated theoretical 
uncertainties. In black are represented the mass radius relations obtained when x = 1 (standard 
GR calculation) while, in red, are represented the mass radius relations obtained without 
self-gravity pressure (namely X = 0). The hatched regions show the theoretical uncertainty 
associated with the fact that the equation of state is in fact unknown. It is clear from the plot 
that this completely dominates the differences between the x = 1 X = 0 situations. Figure 
taken from Ref. [39]. 


where C is an integration constant leading to a curvature term. Now, if x 7^ 1 and 
p = wp, where tc is a constant equation of state parameter, then one obtains 

^2 ^ 1 + 3xw p C 
1 -b ‘iw 3M2 a? ■ 

Using this modified Friedmann equation with tc = 1/3 and ignoring the curvature term 
(which is sub-dominant in presence of radiation as shown by the cosmological data), 
one obtains 


J^2 _ 1 + X Prad 

“ 2 3M2 ■ 


( 1 . 21 ) 


Therefore, the effect of the term proportional to x is to modify the expansion rate of 
the Universe in the radiation dominated era. Or, if one uses the fact that the energy 
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FIG. 3. Light elements abundances calculated when the Friedmann equation is modified 
according to Eq. (1.21). Greens contours are for deuterium abundance, blue ones for helium-4 
and purple ones for lithium-7. The two gray ellipses indicate the region in parameter space 
allowed by observations. Figure taken from Ref. [40]. 


density of radiation is vr^^i^T^/SO, we see that this is also equivalent to changing the 
effective number of relativistic degrees of freedom, namely = g'*(l -|- x)/2. 

Ref. [40] has performed BBN calculations, assuming Eq. (1.21), and computed the 
abundance of deuterium, helium-4 and lithium-7. The isocontours are represented 
in Fig. 3 in the plane {r]io,x)- The parameter gio is defined by rjio = where 

Tj = n^jn^ = pcn/[rn^h?‘2T^C,{J))] ~ 2.73 x [in the last expression, 

TTig ~ 939.6MeV is the baryon (neutron) mass, pcri — 8.099 x 10“^'^GeV^ is the critical 
energy density today, T ~ 2.7255K is the CMB temperature and C(3) — 1.20206]. 
Green contours represent the deuterium abundance {D/H)p, blue contours are helium- 
4 abundance Yp and purple contours are lithium-7 abundance. We see that deuterium 
abundance mainly determines rjiQ while helium-4 abundance gives good constraints on 
the new parameter x- Observations indicate that log {D/H)^ = —4.55 ± 0.04 [41] and 
Yp = 0.24±0.006 [42]. Then, one can identify the region in the space (r/i0)X) which is 
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consistent with those observations. This is indicated in Fig. 3 by the two gray ellipses 
(corresponding to two slightly different assumptions about the abundances inferred 
from the observations). Without entering a detailed discussion, the conclusion is that 
X — 1 is compatible with observations and that the value x = 0 is strongly ruled 
out. Therefore, self-gravity of pressure is, in some sense, confirmed by cosmological 
observations. 

The previous considerations “validate” the mechanism on which inflation is based. 
Inflation thus appears as a well-justified theory. In the next section, we therefore 
describe this theory in more detail and discuss the micro-physics of inflation. 


III. THE MICRO-PHYSICS OF INFLATION OR HOW TO PARAMETRIZE 

INFLATION 

A. Inflation and High Energy Physics 

We have seen in the last section that, in order to have a phase of inflation, we 
need a situation where the fluid dominating the matter content of the Universe has 
a negative pressure. The next question is of course which type of matter can have 
this property. In order to answer this question, let us first remark that inflation is a 
high energy phenomenon by particle physics standards since it is supposed to occur in 
the early Universe. In this situation, the relevant framework to describe matter is not 
fluid mechanics but field theory. And the simplest field, compatible with isotropy and 
homogeneity, is a time dependent scalar field 4>{t) since it has no preferred direction 
and is space-independent. Moreover, in a FLRW Universe, the energy density and 
pressure of a scalar field are given by 

P=y + m, P=\-V{c^). (1.22) 

As a consequence, in a situation where the potential energy dominates over the kinetic 
energy, namely when the field moves slowly or, equivalently, when the potential is flat, 
one obtains a negative pressure and, hence, inflation. The field which drives inflation 
is called the “inflaton”. 
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Let also notice that, when V{(f>) ^ cjP', the equation of state is p ~ —p which, using 
the conservation equation, immediately implies that the energy density, and therefore 
the Hubble parameter H, is almost a constant. The Friedmann equation then leads to 
a scale factor a{t) oc e^L In other words inflation is also a phase of quasi-exponential 
expansion. Moreover, using the expressions established above, one also has 




Pcurv 

Peri 




(1.23) 


and we see that Cljc goes exponentially to zero during inflation. We therefore expect 
to measure a vanishing spatial curvature: this is a first generic prediction of inflation 
and we will see in Sec. V that it is good agreement with the most recent cosmological 
observations. 

As mentioned before, inflation is a high energy phenomenon and, therefore, a 
concrete implementation necessarily rests on high energy physics. In the modern view, 
the micro-physics of inflation should therefore be described by an effective field theory 
characterized by a cutoff A. If the gravitational sector is described by GR, which itself 
is viewed as an effective theory with a cutoff at the Planck scale, then A < Mpi. On 
the other hand, we know that the Hubble parameter during the part of the inflationary 
phase we have observationally access to can be as large as lO^^GeV and this suggests 
A > lO^^GeV. Glearly, at those energy scales, particle physics remains speculative and 
this is the reason why there is currently a plethora of different inflationary scenarios. 
A priori, without any further theoretical guidance, the effective action can therefore 
be written as 


c 

AfpjAg -|— ^ + • • • 

i i 

• • • 5 071; T) -!-••• 



(1.24) 


In the above equation, the first line represents the effective Lagrangian for gravity 
(Ag is the cosmological constant). In practice, we will mainly work with the Einstein- 
Hilbert term only. The second line represents the contribution of matter. We assumed 
that several scalar field are present (a priori, there is no reason to assume that only one 
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field plays a role). The two first terms are the canonical Lagrangian while Oi represents 
a higher order operator of dimension n, > 4, the amplitude of which is determined by 
the coefficient di. Those corrections can modify the potential but also the (standard) 
kinetic term [43]. The last term encodes the interaction between the inflaton fields and 
the rest of the world, i.e. the gauge fields and the fermions T. The dots stand for 
the rest of the terms such as the Lagrangians of of T, the corresponding higher 
order operators etc ... . Notice that the above description is not completely general. 
For instance, suppose that the action of the inflation field is of the Dirac-Born-Infeld 
(DBI) type [44], namely 


S 




R - T{cj>) 



+ r((/))-y(0) , 


where X = —l/2g^^'^d^(l)du4>- An expansion in X gives 

s=y‘d"xv=^(^i?-x+i/+^+•••), 


(1.25) 


(1.26) 


and we see that the higher order terms are not suppressed by a fixed cutoff A but 
by In this case, in some sense, the cutoff has become field dependent. As 

a consequence, the canonical Lagrangian X — 17 is not necessarily always the first 
term of the series and it makes sense to also consider more complicated cases, even at 
“leading order”. 

Another, but related, question is whether the higher order operators can be 
neglected during inflation. Firstly, it is necessary that the field excursion Xcj) be small 
in comparison with the cutoff scale, i.e. Xcj) < A. Whether this is the case or not 
depends on the model. Second, the tree level potential V can receive corrections that 
can be difficult to control. For instance, if there is a mass term, then typically the 
mass m becomes 


m? ^ m? + gM"^ In > (1-27) 

where g is a renormalization scale, M the mass of a heavy field and g the coupling 
between (p and the heavy field. If M > A then one has m > H since we have 
A > H. This means that the potential is no longer flat enough to support inflation, 
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an embarrassing problem indeed! Ways out consist in assuming that the coupling g is 
small or, more convincingly, that symmetries forbid this type of corrections. 

Finally, let us say a few words about the interaction term. Usually, it is considered 
to be negligible during the slow-roll phase. If this is not the case, it leads to warm 
inflation [45—47]. Even if it does not play a role during the accelerated phase, the 
interaction term is of fundamental importance for inflation since it is responsible for 
the reheating stage, that is to say it explains how inflation is smoothly connected to 
the standard hot Big Bang epoch. 

We see that, using theoretical considerations only, it is difficult to restrict the 
Lagrangian of inflation to a simple form. But, in fact, the point is that the CMB 
Planck data can do the job and can constrain the Lagrangian (1.24). For instance, we 
will see in the following that the perturbations are adiabatic (at least for the moment; 
this could of course very well change when more accurate data are collected) and this 
supports the idea that only one scalar field is at play during inflation. Moreover, we 
will also show that Non-Gaussianities have been measured to be compatible with zero 
and this supports the fact that the kinetic term must be standard. We are therefore 
led to consider that inflation is described by the simplest scenario, namely single-field 
slow-roll with a standard kinetic term. It is important to emphasize that we are pushed 
to this class of models, which is clearly easier to analyze than Eq. (1.24), not because 
we want to simplify the scenario but because this is what the CMB data suggest. In 
this framework, the inflationary Lagrangian can be written as 

^ - V{4>) + £int(</>, ^). (1.28) 

In the following, we will ignore the interaction term during the accelerated phase and 
will consider its effect only at the end of inflation (the “reheating” phase). We see that 
we are left with a model that contains only one arbitrary function, the potential V{(j)). 
Therefore, what remains to be done in order to completely characterize inflation is to 
constrain this a priori arbitrary function with cosmological data. This line of research 
has played a dominant role in the recent years. 

Let us now describe the slow-roll formalism which is used in practice to derive 
the inflationary predictions of the models mentioned above. As already remarked 
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previously, one can distinguish two different phases of evolution: the slow-roll phase 
and the reheating phase. In principle, once V{(l)) and '!') are known, the 

model is completely specified. In practice, however, one proceeds in a slightly different 
way. The function V{4>) is considered to be relevant for a limited range of field values 
only, corresponding to our observable window. Then, the evolution of the system is 
controlled by the Friedmann and Klein-Gordon equations, namely 






+ VicP) , 


(1.29) 


+ 3H^ + Vfj) — 0 , 


(1.30) 


where we remind that H = a/a denotes the Hubble parameter and where a subscript (p 
means a derivative with respect to the inflaton field. It is also interesting to introduce 
the Hubble flow functions defined by [48, 49] 


Cn+l 


din \en 
diV 


n > 0, 


(1.31) 


where cq = starts the hierarchy and N = ln(a/aini) is the number of e-folds. 

These parameters provide useful information about the inflationary dynamics. For 
instance, the first slow-roll parameter can be expressed as 


H ^ a 


(1.32) 


and, therefore, inflation (d > 0) occurs if ei < 1. In fact, since the parameters e„ are 
defined in terms of H and since H is determined once F((/>) is known, see Eqs. (1.29) 
and (1.30), it follows that one can also express them in terms of the potential. For 
instance, ei is given by 

ei = -. (1.33) 

2 cP‘^l2 + V{ct>) 

In fact, it is not sufficient to have ei < 1 but one also needs ei <C 1. Indeed, from 
the above expression, we see that this corresponds to a situation where (/>^/2 <C H(^) 
or, in other words, to a situation where the potential is very flat since the field must 
roll very slowly. We just recover the case considered in the previous section. In this 
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situation, referred to as the slow-roll approximation, one has in fact e^, <C 1 for any n. 
If this is the case, then the Hubble flow functions can be expressed as [50] 


— 


62 ^ 2Mp 


e2e3 ~ 2Mp 


V 


y2 


V 


0^00 fV4> 
V (v 


2 / T/ N 4 


(1.34) 

(1.35) 

(1.36) 


The slow-roll approximation allows us to simplify the equations of motion and to 
analytically integrate the inflaton trajectory. Indeed, Eqs. (1-29) and (1.30), which 
control the evolution of 4>, can be rewritten as 


1 + 


^ Mp2(3-ei)’ 

6 - 2ei y dlV ” dy> 


(1.37) 

(1.38) 


As a consequence, in the slow-roll approximation, one has ~ H/(3M|j) and 
(^(f/dN ~ —MpjdlnH/di;!), from which one obtains 


N - iVini = 


v{x) 

Mpi Vxix) 


(1.39) 


^ini being the initial vacuum expectation value of the field. It is clear from the above 
considerations that the inflaton dynamics is entirely determined once the potential 
H(0) has been specified. Since, in addition, the function V{4>) allows us to make 
the connection with high energy physics, it appears as a natural tool to parametrize 
inflation. 


B. Other parametrizations? 

Recently, other parametrizations of inflation have been considered. The motivation 
of these works was to establish a general framework in order to characterize what the 
generic or typical predictions of cosmic inflation are. In this section, we discuss them 
and show that, in fact, they all boil down to choosing a specific potential. 
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The first alternative parametrization that we discuss is the so-called “horizon-flow 
approach” [51-55]. It has been recently discussed in detail in Ref. [56]. Let us define 
a new set of parameters given by 


'A = (2Mp2) 




d(/>^+i ’ 


(1.40) 


Of course, this new definition does not bring any new information. The new parameters 
can be expressed in terms of the previous ones, for instance = ei — € 2 / 2 , 
^A = ef — 86162/2 -|- 6263 / 2 , etc .... It only shows that, if the 6„,’s are all of the 
same order in slow-roll, the ^A are of increasing order. Then, the simple equation, see 
Eq. (1.31) 


is replaced with 


d6n _ 


d6i 

dA^ 

d 62 

diV 

d^A 

d/V 


£ 162 , 


22a - 26? - 36162, 


(1.41) 


(1.42) 

(1.43) 

(1.44) 


The idea is now to truncate this hierarchy at some order M, i.e. to assume that ^A = 0 
for i > M, maybe motivated by the fact that higher order equations deal with higher 
order slow-roll parameters and are thus, in some sense, negligible. Then, this finite set 
of equations (in practice, the case M = 5 has been considered) is numerically integrated 
many times with different initial conditions [51]. In this way, one obtains different 
values of the slow-roll parameters at Hubble radius crossing and, since the observables 
such as the spectral index Uq or the tensor-to-scalar ratio r can be expressed in terms of 
these parameter (see below), different inflationary predictions. The next step consists 
in searching systematic patterns in these predictions which, as a consequence, would be 
considered as “typical” of inflation. In particular, it has been claimed that the different 
predictions for ns and r obtained in this way cluster around the relation [52, 54] 


ri6 ^ ^(1 - ns), 


(1.45) 





22 


J. Martin: The Observational Status of Cosmic Inflation after Planck 


where rig = r/16. The above equation is then viewed as a generic prediction of 
inflation, obtained without the need to specify a particular potential V{4>)- 

However, the above claim is not correct [56]. Indeed, truncating the hierarchy at 
order M clearly means that one assumes, see Eq. (1.40), 


d0M+2 


(1.46) 


an equation which can be easily integrated (!) and leads to [53] 

H{ct)) = Ho 


M+l 

l+'^Ai 

i=l 


Mp 


(1.47) 


Then, from this expression of the Hubble parameter, one can easily calculate the 
corresponding inflationary potential and one obtains 


V{(j)) = 3M2 ^2(0) - 2M^,i7'(,/>) 

M+l 

1+'^ A 


= ^MlHl 


-1 2 


i=l 


Mp 


M+l 

2M3 Mo ^ iAi 


i=l 


Mp 


2—1 


(1.48) 

(1.49) 


The whole procedure is therefore nothing but a particular choice of a potential 
V{(f) depending on M + 1 parameters, Hj. Moreover, Ref. [56] has shown that the 
“mysterious” coefficient 1/3 in Eq. (1.45) can be easily recovered if one carries out a 
standard slow-roll analysis of the potential (1.48). We conclude that this approach is 
not generic at all and only consists in studying a very particular potential. 

More recently, it has also been argued that, rather than choosing a potential V{(f)), 
it is more generic to choose the equation of state during inflation, see Refs. [57-59]. 
So, in practice, what is done is an educated guess for w{N) = p/p. Notice that, since 


l + u;(iV) = -ei(iV), 


(1.50) 


this is also equivalent to choosing a particular function ei(M), which is the strategy 
followed in Refs. [58, 60, 61]. Concretely, one takes 

where a and /3 are two free and positive parameters and Mend is the number of e-folds 
at the end of inflation. However, again, this choice is in fact a choice of V{(f)). Indeed, 
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the slow-roll trajectory (1.39), dN = —Vd(p/{Mp^V') can be re-written as 

2 


d 




'dJV 


dcj) 


dN 


and, from the exact formula 

^ j_/d^y 
2M|i Vdivy ’ 

one obtains the following system of equations 

= 3Ml[l + w{N)], 
^(lnV) = -3[l + wiN)]. 


(1.52) 


(1.53) 


(1.54) 

(1.55) 


When the above set of equations is solved one obtains 4>{N) and V(N) and, eventually 
eliminating N, the function V{(j)). We conclude that giving w{N) and/or €i{N) is not 
a new generic parametrization but just a particular choice of a potential. In order to 
illustrate this point, let us see how it works in practice for the case of Eq. (1.51). The 
trajectory, given by Eq. (1.54), can be written as 


^ = Cl ± -- (iVend - ^ 

iK/pi ol — 


(1.56) 


where Ci is an integration constant. Eor the potential, the integration of Eq. (1.55) is 
also straightforward and one finds 


lnE = C2 + :^(iVend-iV)'"“, 
1 — a 


(1.57) 


where C 2 is another integration constant. Then, from Eq. (1.56), one arrives at 

r a-2 f 4> ^ \ n 2 /( 2 -a) 


(1.58) 


and, inserting this result in Eq. (1.57), one obtains 


hiE = C 2 + 


3/g 

\ — a 


a — 2 / 4> 

■ VP 


-Cl 


2(l-a)/(2-«) 


(1.59) 


This shows that Eq. (1.51) is, in the slow-roll approximation, completely equivalent to 
the choice 


V{(p) = 


( 1 . 60 ) 
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where <5 and 7 are constants [59]. This potential is almost identical to Logamediate 
inflation, LMI in the terminology of Refs. [15, 16], V{4>) = with x = 

{(j)—4>o)/Mpi, in the case where a = 0. This model was studied in detail in Refs. [15, 16]. 
The only difference is that, for LMI, one has a = 4(1 — 7) implying 7 = 1 when a = 0, 
which is not the case here (i.e. a = 0 but 7 is still free). 

We conclude that all the so-called “alternative” parametrizations of inflation 
considered so far are in fact strictly equivalent to specifying a potential. Claiming 
that it is either new or different or better seems dehnitively far-fetched. In addition, 
discussing inflation in terms of V{4>) has the advantage to make the link with high 
energy physics explicit. For these reasons, we conclude that working in terms ofV{4>) 
and scanning the inflationary landscape by considering all possible models seems to be 
the most efficient method to learn about inflation. 


C. Parametrization of Reheating 


Let us now consider the end of inflation, namely the reheating phase, and how 
one can describe it. When ei = 1, the potential is no longer flat enough to support 
an accelerated phase and inflation stops. Usually, this happens in the vicinity of the 
ground state (concretely, the minimum of the potential). At this time, the inflaton field 
starts oscillating and decaying. Then, these decay products thermalize [62] and the 
radiation dominated epoch of the hot Big Bang phase commences. The micro-physics 
of reheating is described by the term £int(<^, T) in Eq. (1.28). But, in fact, in order 
to parametrize reheating, we do not need to have such a detailed description. Indeed, 
as we will see in the following, the inflationary observational predictions are expressed 
in terms of Cn* = Cn (</>*), where iji* is the value of 4> when the pivot scale kp leaves the 
Hubble radius during inflation (the pivot scale is conveniently chosen in the middle of 
the observable window). Since, in the slow-roll approximation, we know the trajectory 
(j) = <f{N), we just need to determine W such that = (/>(W)- This can be done as 
follows. The physical pivot scale during inflation is given by 


kp 

a(iV) 


kp 

®now Oreh Oend ®(A^) 


kp Unow Oreh 
^now ^reh ^end 


(1.61) 
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where Oend denotes the scale factor at the end of inflation and Oreh the scale factor at 
the end of reheating. In the above expression, fep/onow is known and, concretely, we 
take fep/onow = 0.05Mpc“^. The quantity anow/oreh is also known since it only involves 
the standard thermal history of the Universe. On the other hand, the ratio flreh/oend 
depends on what happens during reheating and this is precisely the reason why the 
inflationary predictions are sensitive to this phase of evolution. To go further, we write 
the above equation at the time N = N^. Since, by dehnition, /Cp/a(Ai*) = 

Eq. (1.61) becomes 

i- / V{N:t) _ kp Q-now Qreh 


H{N,) = 


Mp 


(1.62) 


, 3-ei{N,) ttnow ^reh ^end 
the first expression being just the Friedmann equation, see Eq. (1.37). We see that this 
is a transcendental equation for W which, therefore, needs to be solved numerically. 
We also see that it depends on the potential U((/>) and, hence, on the model under 
consideration. Finally, in order to solve this equation, one needs to estimate the 
quantity areh/^end- Let p and p be the total energy density and pressure during 
reheating. Notice that one can have several fluids, possibly interacting which each 
others. The treatment presented here is therefore completely general. Conservation of 
total energy density (we emphasize again that it is not necessary to assume that the 
energy density of each fluid is separately conserved) implies that 


fiV 


p (N) = Pend exp <1 -3 / [1 + Wreh {n)] dn 

J Nend 


(1.63) 


where rcreh = p/p is the “instantaneous” equation of state during reheating. Then, let 
us define the mean equation of state parameter, Wreh, by 


1 


r-^reh 


^reh — 


AN 


Mnd 


Wreh(n)dn, 


(1.64) 


where AN = Nj-eh ~ -^end is the total number of e-folds during reheating. It follows 
that 


and, therefore. 


Preh — PendS 


^AN _ ®reh _ ^ Preh 


-3(1-1-10reh)AW 


®end 


VPend 


-l/(3-|-3?«reh) 


(1.65) 


( 1 . 66 ) 
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As a consequence, the ratio areh/oend depends on two quantities only: the energy 
density at the end of reheating, preh) and the mean equation of state during reheating, 
uJreh- Once a model of inflation is known, pend can be calculated so this is not a 
new quantity (but, again, it introduces an additional dependence on the inflationary 
potential). Inserting Eq. (1.66) into the above expression (1.62) leads to 


HiN,) = 


1 


V{N,) 


n u \ -l/(3+3'U)reh) 
Pend 


(1.67) 


, 3-ei(A^*) ^now 0!-reh 
The above formula still contains Ureh? a quantity that we would like to eliminate from 
the final expression. For this purpose, we write Onow/areh as Onow/oeq x Oeq/oreh, where 
Oeq is the scale factor at matter-radiation equality. Then, we use the fact that, during 
the radiation dominated era, a oc to write 


HiN,) = 


1 


ViN, 




Mp 


3 - ei(iV*) 


/ V Pend / 


aMnd-M 


Peq 

l/4-l/(3-h3^,eh) 


aMnd-M 


( 1 . 68 ) 

(1.69) 


1 /2 

_ ^now Mpi Pgnci Mpi / Preh \ 

Q.now Q-eq pgq^ -^Pl V Pend / 

Except the quantities that are known from standard cosmology (since they only 
depend on post-inflationary physics), such as Onow/aeq x Mpi/plq^, we see that 
this equation singles out the following combination (by definition, the “reheating” 
parameter) [19, 20, 25, 26] 


1/4 

R = ^end p 

^ -O-rad? 


with 


p _ / Preh 

^rad — 1 

V Pend 


Mpi 

-l/4-Hl/(3-H3«J,eh) 


/ Preh \ 
\ Pend / 


(l-3u)reh)/(12 + 12tdreh) 


(1.70) 


(1.71) 


Notice that we have a term Pgnd/-^pi ™ ^9- (1-69). It is introduced because it 
produces a term proportional to the square root of the potential at the end of inflation 
and combines nicely with the vW on the left hand side of Eq. (1.69). The arguments 
presented above can be easily generalized to take into account a change of relativistic 
degrees of freedom between the reheating epoch and today, see Ref. [15]. 

The reheating parameter encodes what can be learned about reheating from the 
CMB. In sec. V, we will see that the Planck data already put constraints on its value. 
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IV. INFLATIONARY PERTURBATIONS 


In this section, we review the theory of inflationary perturbations [63-65] . This part 
of the inflationary scenario is very important because it allows us to use astrophysical 
data to put constraints on cosmic inflation. In the following, we pay special attention 
to the question of how one can calculate the correlation functions of the perturbations 
and to the concept of adiabatic and isocurvature perturbations. As will be seen in 
Sec. V, these quantities carry useful information about the type of inflationary model 
that is realized in Nature. This section can therefore be viewed as a preparation to 
Sec. V in the sense that we discuss in some detail the meaning of the quantities that 
have been measured recently by the Planck experiment. 

To describe CMB anisotropies and large scale structures, one must go beyond 
the cosmological principle. This is a priori a technically difficult task but since 
the inhomogeneities are small in the early Universe, one can use a perturbative 
approximation which, obviously, greatly simplifies the problem. Then, the idea is 
to write the metric tensor as + ''' > where (rj) 

represents the metric tensor of the FLRW Universe and where <C {g)■ 

In fact, 6g^^{r],x) can be expressed in terms of three types of perturbations, scalar, 
vector and tensor. In the context of inflation, only scalar and tensor are important. 
Scalar perturbations are directly coupled to the perturbed stress-energy tensor while 
tensor fluctuations are independent of 5Tfj_y and, in fact, are nothing but gravity waves. 
The equations of motion of each type of fluctuations are given by the perturbed Einstein 
equations, namely = k5T^j_^. 

In order to calculate the behavior of the fluctuations, we also need to specify the 
initial conditions. This is done by postulating that the perturbations are of quantum- 
mechanical origin and that, initially, their quantum state is the vacuum. This is 
possible because, at the beginning of inflation, the physical wavelengths of the Fourier 
modes of the perturbations are smaller than the Hubble radius. This means that, 
initially, space-time curvature is not felt and that, as a consequence, a well-motivated 
vacuum state can be defined. 
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A. Inflationary two-point Correlation Functions 


Once the equations of motion have been derived and the initial conditions specihed, 
one can determine all the statistical properties of the fluctuations, in particular 
their two-point correlation functions or, in Fourier space, power spectra. The scalar 
perturbations are curvature perturbations defined by ({r], £c) = <f> + + <l>)/(3 + 

3w), with w = p/p the equation of state during inflation and being the Bardeen 
potential [64] (not to be confused with the scalar field 4>). As usual in a linear theory, 
it is convenient to work in Fourier space and, therefore, we write 




(1.72) 


As explained before, in the framework of the theory of cosmological perturbations of 
quantum-mechanical origin, the source of the perturbations is the unavoidable zero- 
point vacuum fluctuations. As a consequence, C{r],x) must in fact be viewed as a 
quantum operator and can be expressed as 


Civ,x) 


d^k 

(27r)3/2 




(1.73) 


where and are respectively the annihilation and creation operators satisfying 
[ofc, Op] = 5^^\k—p). The quantum state of the perturbations is the vacuum |0) which 
is, by definition, annihilated by the operator Ofc, namely OfejO) = 0. The function 
gk{p) is the mode function and the Fourier transform of (/{r],x) is given by Cfc(^) = 
akSkiv) + ai^gliri). This last equation leads to (OjCfciCfcalO) = (^i + ^ 2 ). 

From the previous considerations, it follows that the two-point correlation function is 
given by 


where V,^{k) is, by dehnition, the power spectrum of scalar perturbations. An exact 
calculation of this power spectrum is rarely available but a perturbative expansion into 
the slow-roll parameters (since they are small parameters) can be done and results in 


n(fc) 

'PcoiK) 


= a, 


(S) 


+ af^ In 




(1.75) 
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where, as already mentioned, kp is the pivot scale and the overall amplitude can be 
written as 


^0 = 


(1.76) 


Svr^ei^Mpi ’ 

a star meaning that a quantity is evaluated at the time at which the pivot scale 
crossed out the Hubble radius during inflation. We see that the amplitude of the 
power spectrum depends on but also on the hrst slow-roll parameter, ei*. The 
coefficients af'^ can be expressed in terms of the Hubble flow functions. For scalar 
perturbations, at second order in the slow-roll approximation, one gets [48, 49, 66-74] 


ds) _ 


= 1 - 2 (C + 1) ei* - Ce2* + { 2C^ + 2C + — - 5 ] e 
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— —2ei* — € 2 * T 2(2C7 -|- l)e4 T (2C* — l)£i*C2* T ^*62=^ — C^2*^3* T 


(1.77) 

,(1.78) 

(1.79) 

(1.80) 


0-2 ~ “I" 2ei*e2* -|- ^ 2 * ~ ^2*63* + ■ ■ ■ ) 

af^ = 0{el), 

where C' = 7 B + ln2 — 2Ri —0.7296, 7 b being the Euler constant. 

For tensor fluctuations, the approach is exactly similar to what we have just 
described. In particular, the tensor power spectrum Vh can be written in the same 
way as Eq. (1.75) but with a global amplitude now given by 

2H^ 


Vho = 


7r2Mp2 


(1.81) 


This time, the amplitude only depends on the Hubble parameter during inflation. 

Moreover, the coefficients have a similar structure and can be written as 

^2 


dT) _ 


= l-2(C + l)ei. 


2C2 + 2C + y - 5 ) el 


+ ( —C^ — 2C -|- y — 2 ) ei*e2* -|- 


^ — —2ei* -|- 2(2C* -|- l)el — 2{C -\- l)ei*e2* 
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— 4:e4 ~ 2ei*e2* -|- 
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= 0{el) 


(1.82) 

(1.83) 

(1.84) 

(1.85) 
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The coefficients in front of the In k term are related to the spectral indices and, at first 
order in the slow-roll parameters (we will discuss them in more detail in Sec. V), they 
can be expressed as 

ns = l- 2 ei-e 2 , tt-t =-2 ei, (1.86) 

where the first expression refers to scalar perturbations while the second is for tensor 
perturbations. Notice that, sometimes, the power spectrum is written as in the 

context of slow-roll inflation, this is clearly not justified as it would amount to keep 
an infinite number of higher order terms while Uq has been evaluated at first order 
only. It is worth stressing that power-law power spectra are predictions of power-law 
inflation only, that is to say the inflationary model for which V{(f)) oc exp(—C4>) [75]. 
From Eqs. (1.76) and (1.81), one can also estimate the relative contribution of tensor 
and scalar amplitudes 

Vh 

r=^ = lQeu. (1.87) 

' C 

which means that, since ei* <C 1, tensor are sub-dominant. This is of course rather 
unfortunate since a direct measurement of gravity wave would directly lead to the 
energy scale during inflation, 

B. Inflationary three-point Correlation Functions 

We have just derived the slow-roll inflationary two-point correlation functions but, 
of course, higher order correlation functions are also interesting and the field of Non- 
Gaussianity has played an important role in the recent years, see Refs. [76—81] for 
original works on this question and Refs. [82-93] for later works. For a complete 
overview of the subject, we refer to the lecture notes by C. Byrnes [35]. Here, in order 
to be able to fully appreciate the relevance of the Planck data on Non-Gaussianities, 
we discuss how the three-point inflationary correlation functions can be calculated in 
the case of single-field slow-roll inflation with a minimal kinetic term. 

For the two-point correlation, we have seen that it is convenient to work in 
Fourier space and to define the power spectrum. In the same way, for the three- 
point correlation function, we can define the bispectrum as a correlator in Fourier 
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space, namely (Cfci(^) Cfe 2 (^) CfesC^))- feet, we will rather calculate the quantity 
{T^kiiv)T^k 2 {'n)'^k 3 i'n)) where TZ = —'1' — /cf)', 'I' = <1> (valid if a scalar field 

dominates the matter content of the Universe) being another Bardeen potential and 
being the gauge invariant scalar field fluctuation [63]. This amounts to a simple 
change of sign of the three-point function (and no change at the power spectrum level, 
namely V(^ = V-ji, because the power spectrum is quadratic in the Fourier amplitudes) 
since TZ = — Concretely one has 


(7^(?7, x) TZ{ri, x) 7^(r/, x)} = 


d^fci 


d^ko 




(2 7r)3/2 J (2 7r)^/^ J (2 7r)^/^ 


X e 


i (ki+k2+k3)-x 


{TZkiir}) TZk2iv) T^kiij])) 

(1.91) 


In the above expression TZk{ri) obviously represents the Fourier transform of the 
curvature (scalar) perturbation TZ{rj,x), namely 
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TZ{r],x) = 


dkTZkirj) e 


—ikx 


(1.92) 


(27r)3/2 

As explained before, in the framework of the theory of cosmological perturbations of 
quantum-mechanical origin, it is an operator and it can be expressed as 

d^fc 


TZ{ri,x) = 


(27r)^/^ 


akfkivy'"''' + 


t 


—ikx 


(1.93) 


leading to (T^fei'^fca) = l/fci+ ^ 2 ) since TZu = aufk + alfc/fc- Here, the 
creation and annihilation operators are the same as those appearing in Eq. (1.73). 
Of course working in terms of TZkiff) instead of Cfc(^) is both harmless and trivial since 
TZk = —Cfc and /kill) = ~9k{iiy- We do it since many papers on Non-Gaussianities use 
this variable. 

At this stage, it may be useful to say a few words about conventions. In this article, 
we are using Fourier transforms as defined in Eq. (1.92). Another convention, often 


Indeed, the space time component of the perturbed Einstein equation reads 


-('H'h -I- 'h') = k(p -I- p)di 


„(gi) 


(l.g 


(1.89) 


where, for a scalar field, </>'. As a consequence 

1 , Ka^ , , 

cf)' ■ 

Using this last expression in the definition of Z and the Friedmann equation = Ka^pfS, one 
obtains 

C = + = (1.90) 

namely the equation mentioned in the text. 
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used in the literature on Non-Gaussianities, is 

nv, x) = j dk Uk e-*''-", (1.94) 

so that iZk = This implies that {1ZkfT^k2) = (fci + fc 2 ). 

Notice that the two-point correlation function is sometimes defined as (JZkfT^k 2 ) = 
(27r)^P7^(/i;i)5(^) (fci-|-/c 2 ) which leads to the identification Pn{ki) = |/fcj^ [the 
quantity PTiiki) should not be confused with Vniki) = P/(27r^)]. These 

definitions imply that 1Zk = (27r)^/^ {akfk + which can be rewritten as 

Pk = ttkfk + with CLk = (27r)^/^afe. In particular, since [afc,ap] = <5^^^ {k — p), 

we now have [afc,ap] = (27r)^(5^^) {k — p). Different conventions basically correspond 
to different choices for where the factors 2tt appear in the equations. In principle 
straightforward, it can sometimes be confusing when one tries to check a result in the 
existing literature. 

The bispectrum can be evaluated using the standard rules of quantum field theory. 
It is given by [82, 83] 

{TZk^{r])nk2{v)Pk3iv)) = -i / dro(r) ([7^fcl(r/)7^fc2(^?)'l^fc3(f/),i^int(r)]), 

(1.95) 


where r][^i represents an initial time at the beginning of inflation (in practice we take 
ryini —>■ —oo) and r/g a final time at the end of inflation when all the scales relevant 
to the problem are outside the Hubble radius (in practice we take r]e —^ 0). The 
quantity is the interaction Hamiltonian. It can be obtained from the action of the 
system expanded up to third order in TZ, the action of the system being the Einstein- 
Hilbert action plus that of a scalar field (the inflaton). A now standard calculation 
gives [82-84, 86] 


S3[n] = Ml, / dt d^s 


+ ;- ei 62 7t + ^ {d^n) {d,x) {d\) + ^ (5^7^) {dxf + ^ 


elnn^ + aein {dnf - 2 a 61 {d^U) (da) 




2a 


4a 


5C2 


6TZ 


(1.96) 
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where 5 C2I6TZ denotes the variation of the second order action with respect to IZ, and 
is given by 


8 C 2 

5n 


A + HA-ei d‘^n, 


(1.97) 


and the quantities A and x are defined by 


A = 


TZ = a^eiTZ^ d ^A. 


(1.98) 


2M2i?2 

The term T(bL 2 l81^) introduced in Eq. (1.96) stands for the following complicated 
expression 


(8li2\ 


( 8 C 2 \ 


2 a f 8 C 2 


1 


+ 


2aH 


(d^TZ) {da) + 5'^' [A(5*7^) + {d^TZ) {da)] 


X di 




xim^jn 

H- Vj—{di'^) {djTZ) dmdn 

ri 


•-'t) 


(1.99) 


The terms which involves 8 C 2 /87Z can be removed by a suitable field redehnition of 7Z 
of the following form [82-84, 86]: 

7^2 

TZ^TZn + e2^. ( 1 . 100 ) 

After this redehnition, the perturbed action (1.96) becomes a functional of TZn- In the 
following, in order to avoid too complicated notations, we will still use TZ in place of 
TZn- Then, with the redehnition (1.100), the interaction Hamiltonian can be expressed 
as 


Hint{T]) = -Afpi / d^x 


aefTZTZ'^ AaeiTZ {dTZ)^ - 2 ei (9*7^) {da) 


+ ^ ei e '2 (a*7^) {da) {d\) + ^ {d^TZ) {dxf ■ (ITOl) 


where we remind that a prime means a derivative with respect to conformal time. The 
hrst three terms are second order in the slow-roll parameters while the three last ones 
are third order. As a consequence, already at this stage, we see that the bispectrum 
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will be a small quantity. Since we now know the interaction Hamiltonian we can insert 
its expression in Eq. (1.95) in order to derive the bispectrum explicitly. One finds that 

6 


{TZk^irje)TZk2{Ve)T^k3iVe)) = 


(2vr) 




fkliVe) fkiiVe) /fe3(r?e)^c(^l! ^2, ^s) 


+ rkAVe)fk{rie)rk3iVe)Gl{kuk2,ks) 


5 ( 3 ) (fci + k 2 + ks), 


( 1 . 102 ) 


where the delta function ensures momentum conservation. Written in this way, the 
correlator is obviously real. In the above expression, the term Q^{ki,k 2 ,k^) with 
C = (1,6) correspond to the six terms in the interaction Hamiltonian (1.101) (the six 
“vertices”), and are explicitly given by [82] 


Qi{ki, k 2 , ks) 
^ 2 (^ 1 , k2, ks) 
Gsiki, k2, ks) 

Giiki, k2, k^) 
Gbiki, k2, ks) 

Geiki, k2, ks) 


rve 


2i / dr ej (/^^ /fc* + two permutations) 


(1.103) 


'mm 


Ve 


—2i j dra^ f^^ fks (ki • ^2 + two permutations) (1.104) 


l?ini 



+hve permutations , 

rVe 

i dr ei e'a fl^ + two permutations) , 


(1.105) 

(1.106) 
(1.107) 


Vini 



+five permutations , 



+ two permutations 



JL\ 

kl kl) 


{k2 • fca) 


(1.108) 

(1.109) 


( 1 . 110 ) 


Actually, an additional seventh term arises due to the field redefinition (1.100), and 
its contribution to the three point correlation function is found to be 


(^fcl (^e) 77^2 (r^e) 77^3 (r/e))*' ^ — 


{2T^f 62 
{2^f/^ 2 


(l/fczl^ l/feaP + two permutations) 


^ki + k2 + k^). 


(1.111) 
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The other terms in Eq. (1.99) do not contribnte because they all contain a derivative 
(time derivative and/or space derivative) and, at the end of inflation, on super Hubble 
scales, ( = —TZ is constant. 

In order to calculate each of the above terms, one obviously needs to know the 
mode function /fc. Since we evaluate the bispectrum at leading order in slow roll, 
it is in fact sufficient to use the de Sitter mode function, namely fk = iH{l + 
ikr])e~^^^/{2Mpi\/k^ei) (which is properly normalized). Moreover, we only need to 
calculate the first three terms and ^ 7 (^ 1 , ^ 2 ) ^ 3 )) the other contributions being of 
higher orders in slow-roll. In order to illustrate how the calculation proceeds, let us 
explain in detail how ^ 2 (^ 1 , ^ 2 ,^ 3 ) can be calculated (this term is easier than the 
others since we do not have to use the derivative of the mode function). Inserting the 
de Sitter mode function into Eq. (1.104), one obtains 

/ _2 

^ 2 (^ 1 , k 2 , ks) = (^i • ^2 + two permutations) 


X 


dr ^ 


^e t'^(1 - ikiT){l - ik2T){l - ik^r), 


( 1 . 112 ) 


= -2i 


{-iHf 


1 


{ki ■ k 2 + two permutations) 


8M|iy/pppf 

Ve 

1 — ikrpT — {kik2 + ^ 2^3 + kiks) + ikik2k^T^ 






(1.113) 


where k^ = ki + k 2 + k^ is the “total” wave-number. This expression is made of four 
integrals that we need to calculate. The hrst and the fourth ones can be integrated by 
parts and the third one can be directly performed. This leads to 

{-iHf 1 


Q2{ki,k2,kf = —2i 


8Mp3 y/pppf 


(ki ■ k 2 -|- two permutations) 


-1 




Ve 


Vini 


rn 


Ve gifcpT 


- {kik2 + k2k^ -I- kikf) — 


-|- i/Cp J -dr — ik^ 

e‘"'T 


r,ik„ 


T 

Ve 


rVe 

JVini 


Ve gifcpr 


dr 


ikrr 


+ ikik2k‘i 


Vin 


pgifcp? 

ikr^ 


+ 


gifcpi 


Ve 


Vini 


, (1.114) 


and we see that the second integral exactly cancels the term arising from the integration 
by parts of the first integral. In principle, at this stage, it is sufficient to take rji^i = —00 
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in the above expression in order to get the final result. But, obviously, the result would 
be ill-defined. So what is done is to slightly rotate the integration path in the complex 
plane and replace r/ini with —oo(l — i6) where 5 is a small parameter. This produces 
a term which, in fact, kills all terms proportional to It is worth 

noticing that this should not be viewed as an arbitrary technical trick but as the 
standard method to properly identify the correct vacuum state [94]. As a result, one 
obtains the following expression 


Q2{ki,k2, fcs) = —2i 


{-iHf 




e?* (^1 • ^2 + two permutations) 




+ fesfes + felfcs) 

rje ik. 


+ikik2k3 




ikr£ Tje 


ikrr 


+ 


^ikr^Tje 


kl 


(1.115) 


Then, the final step is to take r/e —)• 0. Clearly, there is a problem with the first term 
and, therefore, in the following expressions, we will keep r/e unspecified. For the other 
terms, the above expression simplifies and one is led to 


G2{ki,k2, fcs) = —2i 


{-iHf 


1 


8Mp3y^pppf^2 


e?* (^1 • ^2 + two permutations) 


-1 


— ^ ikik2 + k2k3 + kiks) + ^kik2k3 


Ve 


k, 


k^ 

T 


(1.116) 


This completes the calculation of ^ 2 (^ 1 , ^ 2 ) ^a)- Now, we insert the above re¬ 
sult into Eq. (1.102) in order to determine the contribution of ^ 2(^11 ^ 2 ) ^ 3 ) to 
{T^kiiVe)T^k 2 i'ne)T^k 3 {'ne))■ To perform this calculation, we need fkiVe), which we 
take to be iH/[8Mp^^/€^{rje)kfkf^] since the limit r/e —)• 0 does not cause any problem 
in that case. As a result, one finds that 


(7^fc,(r/e)7^fc,(r/e)7^fe3(r/e))(2) = (fc^ + + ^ 3 ) 


{iHf 


_ 8M|i 


G2{ki,k2, ks) + 


{-iHf 


8M|iY^ef(r/^)A!pp: 


,g*^{k3,k2,k3) 


(1.117) 
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which, combined with Eq. (1.116), leads to 


^ ■ ^2 + two permutations) -^(-^) i — 

8 Mp 3 H i I 


^ikr^ fje 


% z 

+ — {kik2 + k2k3 + kiks) + -r^kik2kz 

rtrp rC^ 

% % 

- T^{kik2 + k2k3 + kiks) - ■^/ciA:2^3 

/Ct f^rj. 




-1 

Ve 




(1.118) 


This expression can be simplified further and one obtains the following formula 


(7^fc,(r?e)^A,.(r?e)7^fc3(%))^') = ^^|^MpV(')(fci + fc2 + fc3) 




QAHm^,et\\/\r^,)klklkl 

sin(A:p7?e) , 2 


(fei • k 2 + two permutations) 


(/ciA :2 + k 2 kz + /cifcs) + 7^^1/22/23 


/2.p77e /2 t: 




(1.119) 


We see that the limit ? 7 e —t 0 is now well defined and can be taken. The term in 
Q 2 {ki, k 2 , k^) was singular but, combined with its complex conjugate in the correlator, 
the limit has become regular. Therefore, the appearance of a singular limit was just 
a temporary technical problem and, in the expression of the physical quantity, the 
problematic term has disappeared. The final expression reads 


("/^fcl (%)'/^fc 2 (^e) ^fcs (??e))^ ^ ~ /rv \Q 




(27r)9/2 IGM^ei {kik2k^Y 

-/2p + (/21/22 + /22/23 + /21/23) + Ty/21/22/23 

rCrp rC^ 


{ki ■ k 2 + two permutations) 


5 ( 3 ) {ki + k2 + k^ 


( 1 . 120 ) 


As expected the amplitude is controlled by the Hubble parameter (to the power four 
while the amplitude of the power spectrum was quadratic in H) and the (first) slow-roll 
parameter. We also see that the scale dependence is quite complicated. 

The calculation proceeds exactly the same way for the first and third terms. 
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Explicitly, one obtains 




1 




/urp J /trp 


( 27 r) 9/2 IGM^ei {kik 2 hf 

^ klkj 

krj, J krj, V krj, J k rj. 


5(3) (fci + k2 + fcs) , 
( 1 . 121 ) 




1 


(7^fc,(r/e)7^fc2(%)^fc3(%))^^^ (27r)9/2 16MpVi (A^ifcaA^s)^ 


(fcl • k2)-r^ ( 2 + 

/X/r-T- 


fci + A:2 
/c™ 


+ (A^2-fc3)^ (2+^^- 

/Vrr 


A-2 / 

+ (^1 ■ ^ 3 )^ 1^2 + 
5 ( 3 ) (fc^ + fc 2 + fcg) . 


ki + fc3 


( 1 . 122 ) 


l/rp Y »l.rp 

Finally, the seventh term given by Eq. (1.111) can be re-written in terms of the two- 
point correlation function 


(■^fcl(%)’^fe2(%)'^fc3(^e))(^^ = l^^n/o 27r^e2 


1 


(27r)9/2 {kik2k^Y 

k\Vn{k2)Vn{k^) + klVn{ki)Vn{k^) + klVn{ki)Vn{k2) <^(3) {k^ + fca + fcs), 

(1.123) 


where, in order to evaluate the last term, we have made use of the definition introduced 
before: Vnik) = A:3|/;^,|2/(27r^), see Eq.(1.74). 

We have now completed the calculation of the three-point correlation function in 
Fourier space. We notice that, as already mentioned above, the dependence in ki, k 2 , 
ks is rather non trivial. In order to emphasize this point, it is interesting to recalculate 
the three-point correlation function in the following simple setup. Suppose that we 
write the curvature perturbation as 

3 

n{ri,x) =n^{r],x) - ^'R,l{rj,x) + -- - , (1.124) 

where TZ^ denotes a Gaussian quantity, and the factor of 3/5 arises due to the relation 
between the Bardeen potential and the curvature perturbation during the matter 
dominated epoch. The amplitude of the quadratic term is constant and conventionally 
called f^^. Let us notice that this assumption is highly non trivial and that, a priori, 
the coefficient in front of the quadratic term is expected to be a function of space. 
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Postulating that it is a constant enforces a particular scale dependence of the three- 
point correlation function as we are going to see. In Fourier space, the Gaussian part 
is written IZ^ = (27r)“^/^ J dfcT^^e"*^'* and it follows that 

= (2^ / dfc(27r)-3/2y‘ dp7^«7^°_pe-*^■^ (1.125) 

from which we can read the Fourier coefficient of the non-linear curvature perturbation, 
namely 

3 f 

7^fc = - ^(27r)-3/2 j dp7^G (1.126) 

Using this expression, one can now evaluate the bispectrum. One obtains 

/ r 3 r 

= ^(27r)-3/2y dpl7^G 

O ^loc /> 

X n, - -^{2^)-^'^ J dp2 K$, 

X n,-^^2„)-^'\l dpzK,n,-v:. (1127) 

and, therefore, 

-^(27r)-3/2 J dp 3 (i^fc^(??)1^fe2(r/)’7^)l3(?/)i^fc3_P3(?/)) + two permutations -b • • • , 

(1.128) 

where the dots denote the higher order terms. Since the three point correlation function 
vanishes for Gaussian statistics, the previous expression reduces to 

(JZuMT^kMT^kM) = -^(27r)-3/2y dp3 

+ two permutations (1.129) 


As expected, the three-point correlation function is proportional to the coefficient 
To proceed, one can evaluate this expression by means of the Wick’s theorem. Then, 
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one obtains 


(7^fc,(r/)7^fc2(r/)7^fc3(r/)) = 
3 .I 0 C 


-(2vr: 


1 - 3/2 


dp3 






+ two permutations + 


(1.130) 


Since the two-point correlation functions are nothing but the power spectrum, the 
above expression takes the following form 


(7^fc,(r/)7^fc2(r/)7^fc3(r/)) = - 


3/r 


5 

\2 - 


(27r: 


1 - 3/2 


ik, + fc,) (^^,(3) (^3 + ^3 _ p3) 

+ (fc, +P3) ,,,3 + _ P3) 


X J dp3 

2 /cf - 2 

'(2 


uS 

^2 




-|- two permutations -|- 


u3 

^2 


(1.131) 


Then, the integral over can be easily performed, thanks to the presence of the Dirac 
delta functions. We see that the first term in the above expression is different from the 
two next ones. Indeed, it leads to a term which can be ignored since, in some 

sense, it is homogeneous and only participates to the background. The two other terms 
yield a (5^3) -^ 1^2 + ks) which ensures momentum conservation. The final expression 
reads 

3 1 

{nkfRkfRk,) = -^ (2 7r)4 (2 7r)-3/2 ^3 ^3 ^3 ^^"^(^1 + ^2 + ^ 3 ) 

X [^1 'Pn{k 2 ) 'Pniks) -I-two permutations] . (1.132) 


We see that the scale dependence of the bispectrum for this simple model does 
not reproduce what we obtained in the case of inflation, see Eqs. (1.120), (1.121), 
(1.122) and (1.123). The inflationary case is clearly much more complicated. In fact, 
Eq. (1.132) has a similar structure as {TZkg{r]e)'R-k 2 i'ne)T^k 3 {Ve))^'^\ see Eq. (1.123). 
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But the three extra terms (r/e )’^^2 (^e)are such that the full slow- 
roll bispectrum differs from Eq. (1.132). 

At this stage, it is worth discussing again our conventions. We have seen below 
Eq. (1.94) that, often in the literature, the two-point correlation function is defined 
as (■^fci'^fca) = (27r)3p7j(fci)(5(3) (fci-k ^ 2 ), where = ( 27 ^) 3 / 27 ^^^ and Pn{ki) = 
l/fciP 7 ^ Pniki). Then, in order to mimic and/or generalize the definition of the two- 
point correlation function, the following definition of the bispectrum B-ji{ki,k 2 ,k 3 ) is 
introduced 

{1Zki'Pk2'Pki) = (2vr)^B7e(fci, /c2, fcs) (fci 4-^2 + ^ 3 ) • (1.133) 

Notice that we could have also used another definition {'JZki'Pk 2 'Pk 3 ) = 
(27r)^B7^(A:i, A: 2 , fes) (fci-|-A :2 + ^ 3)5 which would have resulted in a difference 
by a factor of (27r)®/^ [and, by the way, explains the appearance of such a factor in 
Eq. (1.102)]. Here, we do not follow this route and use the convention (1.133). Then, 
Eq. (1.132) implies that 

Bn{ki,k 2 ,k 3 ) = (l/fcal Vfcgl^ + two permutations) (1.134) 

= -^/nl [Pn{k 2 )Pn{k 3 ) + two permutations] . (1.135) 

It is also frequent to define the bispectrum of the Bardeen potential rather than the 
conserved quantity TZ. Concretely, the definition reads^ 

{^ki^k2^k3) = (27r)^B$(A:i, ^ 2 , ^ 3 ) (^i + ^2 + %) • (1.138) 

Since ( = 54>/3 = —P, we have B^ = —27B'ii/125. However, since fk = —{5/3)uk, we 
also have P-ji = (25/9)P$. As a consequence, from Eq. (1.135), one obtains that 

H$(A;i, k 2 , fca) = 2,f^_^ [P^{k 2 )P^{k 3 ) + two permutations] , (1.139) 

^ As already mentioned, our convention for the Fourier transform is such that 

^iv,x) = J dfc4>fc(?7)e~^'° °’ (1.136) 

and, following the notation that we have already introduced, = {‘2Tv)^^^^k- Moreover, if the 
Bardeen potential quantum operator is written as 

l-(77, x) = J [afcUfc(77)e”'"“ -f 4 u^77)e"“'“ ”] , 

then one has P# = \uk\^ and = k^\uk\‘^/{2P). 


(1.137) 
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which is a formula that often appears in the literature. 

Of course, we can also put the slow-roll bispectrum calculated before under the 
form given by Eq. (1.133). For this purpose, let us write Eqs. (1.121), (1.120), (1.122) 
and (1.123) as 

(7^fcl(%)7^fe,(r?e)7^fc3 (%))(') = + ^2 + fc3), (1.140) 

where the concrete expression of the can be read off from those equations. Then, 
the bispectrum for single-field slow-roll models can be written as 




1 

Z^O Z-aO 1^0 

^2 ^3 


E 

i=l,2,3,7 


(1.141) 


The previous result can also be used to define an effective, scale dependent, /j.,l 
parameter. If we equate the full bispectrum X]i=i 2 3 7 (^fci (^e)'l^fe 2 (^e)'l^fc3(7/e))^*^ to 
the expression of Eq. (1.132), one obtains 

1 = 1 , 2 , 3,7 

X [fc?’P 7 e(^ 2 ) T’ 7 e(A: 3 )-|-two permutations] (1.142) 


If, for instance, we evaluate this quantity for ki = —^2 and a vanishing (so that 
ki + k 2 + k^ is also zero which is mandatory given the presence of the Dirac function 
in the above expressions), then the expressions of simplify such that one obtains 


E 

1=1,2, 3,7 


/I 3 62 

IGM^iCi V2 ^ 2 ^ ^ ^ 


fP^k^ 


(1.143) 


which, using Eq. (1.86), can be written as [82, 95] 


where “sq” means “squeezed” and refers to the fact that we have taken the particular 
configuration ki = —k 2 and a vanishing k^. Notice that, since we calculate a three- 
point correlation function, the sign of fj^^ is non trivial. The sign that we have 
obtained results from the choice made in Eq. (1.124) and from the fact that we evaluate 
the correlator of TZ. Finally, very roughly speaking (see Sec. V for a more complete 
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discussion) the present status of the art is such that one can detect Non-Gaussianities 
if > 5. For slow-roll models, since Ug — 0.96, one obtains — —1.6 x 10“^, 
a number that is therefore undetectable. This conclusion is in fact valid for any 
configuration one may choose. Let us also mention that other consistency relations for 
Non-Gaussianities have recently been studied in Refs. [90, 92, 93]. 

Clearly, a detection of a non-vanishing three-point correlation function, given 
present day technology, would immediately rule out single field slow-roll models with a 
standard kinetic term. It is therefore quite remarkable that Non-Gaussianity has not 
been detected so far. Let us also stress that the opposite statement is not true. The 
fact that we do not see Non-Gaussianities does not imply that the more complicated 
models of inflation are necessarily ruled out even if some of them do predict large 
Non-Gaussianities. For the calculation of the three-point correlation functions of these 
more complicated models, we again refer to Ref. [35]. 

C. Inflationary four-point Correlation Functions 

Obviously, the next step is to calculate the four-point correlation function or 
trispectrum [96-98]. Of course, when we consider higher order correlation functions, 
the calculations become more and more complicated. In the previous sub-section, 
we calculated the action at third order in the perturbations in order to derive the 
inflationary three-point correlation function. In order to calculate the four-point 
correlation function, one therefore needs to evaluate the perturbed action at fourth 
order. In order to get an idea of how involved it can be, let us consider again Eq. (1.124) 
but expanded up to third order 

n{rj,x) = n^{r],x) - '^^TZl{r],x) + x) ■ ■ ■ , (1.145) 

thus introducing the parameter Here, we write in order to avoid cumbersome 
notations but it should be clear that (and this will be the case in the rest 

of this section). The cube of the curvature perturbation can be expressed as 

7^^(r/,®) = j dk{27r)-^ j dpdq7^fc_p_q7^p7^qe"*'“■“’, 


(1.146) 
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which allows us to identify the Fourier transform of the cube of the curvature 
perturbation [as we identified the Fourier transform of the square of the curvature 
perturbation in Eq. (1.126)]. Then, the four-point correlator takes the form 


(JZki'R-k2'T^kfTlki) = 

3 /nl 


-(2vr)- 


n,- 5 




dP2W» + 


(2»)-3/'^ydp4 7J» + ?|i(2») 



|dpldgl7^G 7^G 


dp2dq2ng_^^_^n^^n°^ 


dp3dq3TZl^_p^_^n%n% 


dp4dq4^g,_p,_g^^®4^°4 


(1.147) 


Expanding this expression, one arrives at the following formula 


{nk^TZk^nk^TZk,) = {nin^n.n,) 

- ^^( 27 r)- 3 / 2 (^y'dp 4 (^fei^fc2^fc3^p4^fe4-P4) + permutations^ 

+ dpi j dp2 ( 7 ^° + five permutations^ 

+ |5.j2vr)-3(/ dp4 1 

+ three permutations^ + • • • , (1.148) 

where the dots denote higher order terms. The first term in the above expansion is 
non-vanishing but can be expressed as the square of two-point correlation functions 
and will be ignored in the following. The second term is zero since it involves five- 
point correlation functions of Gaussian quantities. The two last terms are the terms 
of interest. We see that they are given in terms of a six-point correlation function, a 
quantity which is not zero for Gaussian quantities. These terms can be evaluated by 
means of the Wick’s theorem and lead to the sum of fifteen terms, each of them being 
made of the product of three two-point correlation functions. Eor the term proportional 
to f '^^, among the fifteen only eight of them actually contribute. An example of a term 
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contributing is given by 

/dp,dp, <w» {n.-pM.) 

^ /dpidp2i^^7’K(pi)d*^tpi +P2)^^^'PK(|fcl -pil)d*^'{^i -Pi +^ 3 ) 

X i^^'PK(|fc2-P2l)d^^t^2-P2 + fc4) 

O J fcg 

X d(3)(fc2 + Pi + ki) 

= |fc^+fc3|3fc3fc3 ^^(l^l + fc3|)iP7^(fc3)^7^(A:4)^("^(fcl +^2 + ^3 + ^ 4 )- (1.149) 

In fact among the eight terms mentioned above, four are identical to the one we have 
just calculated and the remaining four are all given by Eq. (1.149), but with |A:i + ^ 3 ! 
replaced with |fci + ki\. On the other hand, an example of a non-contributing term is 

I dp.dp2 {ni,n,-p,) i'K-i.n.-p.) {nM.) ■ (‘-iso) 

We see that the first two-point correlation function appearing in the above integral 
will lead to a term proportional to + ki — pi) = 5^^\ki), which, in some sense, 

is homogeneous. This explains why the term in Eq. (1.150) can be ignored. 

Let us now come back to Eq. (1.148) and consider the term proportional to 
Using again Wick’s theorem, this term can be expressed as the sum of fifteen 
terms made of the product of three two-point correlation functions. Among these 
fifteen terms, only six participate to the final expression (and they all give the same 
contribution). One example is 

J dp4dq4 

= J dpidq4^^^VTz{h)S^^\ki + k4-pi-qi)^^I^V'ji{k2)S^^\k2+Pi) 

X ^2^3 'P'R.{k'i)5^'^\k2, + 94 ) 

= ^ o , 3 , 3 , 3 '^ 7 ^(fel)'P 7 ^(fe 2 )L^ 7 ^(fe 3 )<^^^^(fcl + k 2 + ks + k^). (1.151) 

o /C^ /C 2 
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Putting everything together, one obtains the following expression 


9 


4 X 


(2vr) 


1 


+eleven permutations 


8 |fci + fesP/cg/cl 
9ff: 


rn{\ki + k3\)Vn{k3)Vn{ki) 


+ k2 + ks + ki) + ^ 

25 


6 X 


(2vr)6 1 

o 1.31.3 h3 


^'Pn{ki)'Pn{k 2 )'Pn{k 3 ) + three permutations 

36/^ 


25 


5^^\ki + k2 + ks + k^) (1.152) 

(27r)“^( |/fci+fc 3 p|/fe 3 p|/fc 4 p + eleven permutations )+ ^2 + ^3 + ^ 4 ) 


+ - 


^^^(27r) ^f|/fciPl/fc 3 l^|/fc 3 p + three permutations^ (5^^^(fci + ^2 + ^3 + ^ 4 )- 


25 


(1.153) 


The fact that we have eleven permutation in the first term comes from the fact that 
we had six terms and that each of these terms separates in two groups. At the end, 
this gives twelve terms. Usually, the definition of the trispectrum is given in terms of 
iZk: (see the above discussions about conventions) and reads 

= (27r)^77e(A;i,/c2, ^3) fc4)<5^^^(fcl + ^2 + ^3 + ^ 4 ), (1.154) 


with 


Tn{ki,k 2 , k 3 , k^) = ( |/fci+fe 3 l^l/fe 3 l^l/fe 4 l^ + eleven permutationsj 


+ 


25 


l/fcil^l/fcsl^l/fcsl^ + three permutations ). (1.155) 


One can check that our result (1.153) matches exactly this form provided that 

36/^,, 


25 


(1.156) 


This equation is called the Suyama-Yamaguchi consistency relation [99] (more precisely, 
it is in fact a particular case of ^ 36/^^/25). The above equation indicates that 
the tri-spectrum is expected to be quadratic in the slow-roll parameters and, hence, 
even harder to detect than the three-point correlation function. Of course, it should be 
stressed again that the scale dependence of Eq. (1.155) is not what would emerge from 
an exact calculation starting from the perturbed action at fourth order. In Sec. V, we 
will discuss the constraints put by the Planck experiment on the tri-spectrum. 
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D. Adiabatic and Isocurvature Perturbations 


Another important consequence that follows from the Planck data is that the 
perturbations are adiabatic. Before discussing in more detail in Sec. V how this 
conclusion is reached, we now explain what it means and what it implies for inflation. 

The post inflationary Universe is made of four fluids: photons, neutrinos, baryons 
and cold dark matter (we are ignoring dark energy). In order to calculate the CMB 
anisotropies, one needs to integrate the equations governing the behavior of these 
four fluids. But we also need to specify initial conditions, just after inflation, at the 
onset of the radiation dominated era. Different initial conditions will lead to different 
subsequent evolutions and, therefore, to different CMB patterns. Adiabaticity refers 
to a situation where one has [ 100 ] 

<5cdm = <^b = ^<^7 = (1.157) 

where 6 x = Spx/px is the density contrast (“cdm” stands for cold dark matter, “b” 
for baryons, 7 for photons and ly for neutrinos). It may be surprising that CMB 
data single out particular initial conditions and it is interesting to discuss why the 
conditions (1.157) play an important role. Equally important is the question of what 
they can teach us about inflation: after all, these initial conditions are the results of 
what happened during inflation. As a consequence, they certainly tell us something 
about the type of inflationary expansion that took place in the early Universe. 

Let us start by giving the equations controlling the evolution of the four fluids 
mentioned before. Each fluid is characterized by its density contrast 6 x and by its 
velocity vx- Erom energy conservation, one can derive the following equations 


((5c — ST)' — k^Vc = A(, — k'^Vc = 0, 

(1.158) 

((5b - 3T)' - = A'b - fe^Ub = 0, 

(1.159) 

- 4T)' - - ^k‘^v^ = 0, 

0 0 

(1.160) 

- 4T)' - = 0, 

(1.161) 


where T is the second Bardeen potential already considered before (but, in the present 
context, we no longer necessarily have T = <I>) and where the quantities Ax are defined 
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by the above equations. The space component of the conservation equation gives an 
equation for the velocities. For cold dark matter, one obtains 

v'^ + nvc + ^ = 0, (1.162) 


where is the other Bardeen potential. In the early Universe, baryons and photons 
are tightly coupled. This means that Vh = = Vby- The corresponding equation of 

motion reads 


^b7 + 


R 


1 + R 


T *1? + 


1 5^ 


+ 


4?7b7 R 


A1 + R 3ap^ 1 + i? 


= 0 , 


(1.163) 


where is the viscosity (or anisotropic stress) of the fluid made of baryons and 
photons and R is three quarters of the baryon to photon energy density ratio, namely 
R = 3p\^/{4:p^). Finally, the conservation equation for the neutrinos can be written as 


^ + -tSu H—— k^Vi, — 0, 

4 ap„ 


(1.164) 


where pu is the neutrinos viscosity (notice that the viscosity does not appear in the 
time component of the conservation equations). Since the above formulas contain the 
two Bardeen potentials, they must be supplemented by additional equations governing 
the behavior of and T. These are of course the perturbed Einstein equations. By 
combining the time-time and time-space Einstein equations, one arrives at 

+ = + wx)vx, (1.165) 

where the sum runs over the four species mentioned above, where wx is the equation 
of state parameter of the fluid X and Tf = a'/a. Einally the space space component 
of the Einstein equations (with i j) leads to 


^2 

— ($ - ^r) = (?7b7A:^fb7 + Puk'^Vy) 




Q-pcri 


(1.166) 


where we remind that pcri is the critical energy density. At this stage we have all 
the equations necessary to understand the behavior of the four fluids: we have ten 
quantities (namely four 6x, four vx, 'k and <k) and ten equations, namely Eqs. (1.158), 
(1.159), (1.160), (1.161), (1.162), (1.163), (1.164), (1.165) and (1.166) (the tenth 
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equation is simply Ub = v^). The only thing which remains to be done is to specify the 
initial conditions. Integrating this system of ten equations analytically is not possible 
(even if linear). This has to be done numerically. However, since we are mainly 
interested in the behavior of the system on large scales, the problem gets simplified. 
Indeed, let us introduce the quantity, introduced by Bardeen, Steinhardt and Turner, 
Cggrp defined by [101, 102] 

^ X ^ 

where p = px is the total energy density and (x can be expressed as 


Cx 




3(1 + rex)’ 


(1.168) 


From Eqs. (1.158), (1.159), (1.160) and (1.161), we see that, on large scales (where the 
terms oc k'^vx go to zero), each (x is conserved, namely C'x — 0- Now, we understand 
the particular role of the conditions (1.157). Indeed, they amount to simply choose 


Cedm — Cb — C7 — Cu — Cadia- 


(1.169) 


and, in this case, we have 


Cbst — Cadia ^ ^ , — Cadia j 

X ^ 


(1.170) 


which is a constant. Therefore, for adiabatic initial conditions, the quantity Cbst 
conserved on large scales. Another way to see the same thing is to differentiate Cbst 
( using the expression of 6 p' obtained from energy conservation). Then, one arrives at 
the following equation 


CL, = -^<5Pnad - , (1.171) 

BSl p + P 3 

which shows that, on large scales, the conservation of Cbst controlled by the non- 
adiabatic pressure [here, is the scalar component of the gauge-invariant velocity]. 
This quantity is defined by the following expression 


<5pnad = 6p- c^dp 


(1.172) 
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where 6p, 6p are the total perturbed energy density and pressure, respectively. The 
quantity = p' jp' is the (total) sound velocity. In the case where one has two fluids 
(in order to keep things simple), expressing the perturbed energy density and the 
perturbed pressure explicitly, one arrives at 


(Jpnad = {5pi - Cg/pi) + {Sp2 - cIJP2) 


,( 2 2 \ +Pl) (P2 +P2) c, 

+ Msi “ S 2 ; --*^12 > 


where 5i2 is given by 


^12 = 


bp\ 


bp2 


Pi +Pl P2+ P2 


p + p 


— 3 (Cl - C 2 ) • 


(1.173) 


(1.174) 


and where Cg; = p^p'^- The non-adiabatic pressure contains two contributions. The 
terms 6pi — Cg.hpi originate from intrinsic entropy perturbations (if any) of the fluids 
under consideration while the term proportional to S 12 represents the entropy of 
mixing. Let us summarize: for adiabatic perturbations, Cbst ^ conserved quantity. 
For non adiabatic perturbations, this quantity can evolve even on large scales and this 
evolution is given by Eq. (1.167). 

Let us also remark that one can work in terms of the quantity C defined by [101] 
and already introduced before 


c 


4> + 


3 1 +UJ 


(1.175) 


If one has then 


^BST 


-c 


a. 

-$ 

3ei772 ’ 


(1.176) 


and, in the standard situation, when there is no entropy perturbations, the quantities 
C and Cbst both conserved on super-Hubble scales. Notice that, strictly speaking, 
C stays constant only in absence of shear viscosity. 

Let us now try to understand how the presence or the absence of adiabatic 
perturbations can affect CMB anisotropies. On large scales, the temperature 
fluctuations can be expressed as 

1 . , 

JP ~ 4'^7llss T 4?|iss, 


(1.177) 
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where “Iss” means “last scattering surface” and indicates when the radiation density 
contrast and the Bardeen potential must be evaluated. Since last scatterings occur 
during the matter dominated era, using the time-time component of the perturbed 
Einstein equation, one obtains —2<l>|iss ~ .Rcdm<Jcdm|iss + -Rb<5b|iss where i?cdm = 
PcAra/{Pcdm + Ph) and = Ph/{PcAm + Ph)- It is Conventional to measure the non- 
adiabatic perturbation with respect to photons. Therefore, one introduces the notation 
Sx = Sx'y = 3(Cx — C 7 )- Then, one obtains 

•^llss = “ F-^cdm^cdm “ (1.178) 

5 5 5 


where we recall that i?b and i?cdm are evaluated at last scattering. Notice also 
that, in principle, we do not need a subscript “Iss” for Qx or Sx because they are 
constant (in time) quantities since C,x is conserved. In particular, they should be 
viewed as the value of Cx at the onset of the radiation dominated era, just after 
inflation and, therefore, Sx could also be written as S'™ in order to emphasize 
this point. To calculate the temperature anisotropies, we use Eq. (1.177) and write 
ST/T = (5.y/4|iss + ‘hiss = C 7 + 2‘h|iss- As a result, one obtains 


ST 

~¥ 



2 

5 


77cdmScdm 


2 

5 


RhSb- 


(1.179) 


Finally, during the Radiation Dominated (RD) era, one can write 


Su 

Crd ~ R7C7 T RuCu = C7 T 


(1.180) 


with = p^/{p^ + py) and Ry = py/{p^ + py), these quantities being evaluated during 
the radiation dominated era. Using Eq. (1.180) to obtain an expression of and using 
this expression in Eq. (1.179), it follows that 


ST_ 

~¥ 


r Crd r I^cdm'S'cdm r T 

5 5 5 15 


(1.181) 


which coincides with Eq. (7) of Ref. [103]. The term Crd/^ represents the adiabatic 
contribution. In fact, one can also define an effective isocurvature mode taking into 
account both cold dark matter and baryons entropy fluctuations by dehning 
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such that Eq. (1.181) now reads 

— ~^Crd ~ g-^cdm'S’cdm + (1.183) 

We therefore have two adiabatic modes that are, as will be seen in Sec. V, denoted 
by the Planck collaboration GDI (for the effective cold dark matter) and NDI (for 
neutrinos). In fact, there is a third mode, NVI, related to neutrinos velocity. Since 
the expression of the temperature is modified by the presence of isocurvature modes, 
the temperature multipole moments will also be affected, for concrete and quantitative 
results see for instance Ref. [104]. As a consequence, when compared to the CMB data, 
one can put constraints on their amplitude. 

As will be discussed in Sec. V, so far, CMB measurements are consistent with 
adiabaticity. This gives non trivial information about inflation. Indeed, if non adiabatic 
pertubations were observed it would mean that inflation can not be driven by a single 
scalar field. As for Non-Gaussianities, this would have implied that single-field slow- 
roll inflation with a standard kinetic term were ruled out. This class of models has 
therefore passed another non-trivial test. Of course, this is the situation now and this 
could very well change in the future. In that case, what would be the implications for 
inflation? A natural explanation would be to have multiple field inflation and we now 
explain in detail why using a simple example [103-115]. 

Assume that, instead of having one field, we now have a collection of fields 
that all play a role during inflation. For simplicity, and because we want to be 
explicit, let us consider the case where we have two fields, (/>h and (pi, and where the 
potential is quadratic for each field, without interaction term, namely V = I 4 + Rf = 
m‘^<p‘^j2 + m‘f(pi/2. Then, the equations of motion for the background are given by 



(1.184) 

0h + 3F(^h + mlcph = 0, 

(1.185) 

(pi + 2>H(pi -)- m\(pi = 0, 

(1.186) 


where, as is standard in the literature, we have used notations that make obvious the 
fact that one field is heavy and the other light, meaning that R = mh/rrii > 1 [not to 



Cosmic Microwave Background Radiation: II JPBCosmo 


53 


be confused with the R introduced in Eq. (1.163)]. These equations cannot be solved 
exactly but one can use the slow-roll approximation. The first Hubble flow parameter 
is given by, see Eq. (1.32) 


{i’h + ) (1.187) 

and ei <C 1 implies that <C 1 and K(^|/(277^) <C 1. These two conditions 

are similar to what would be obtained in the single-field case. This means that, as 
usual, the kinetic term can be neglected in the Eriedmann equation. On the other 
hand, the second Hubble flow parameter can be written as 


62 = 2ei -I- 


2 (/>h0h + 


(1.188) 


H </)^ + 

In the single-field case, this relation reduces to 62 = 2ei + 2(j)/(Hep) and the acceleration 
in the Klein-Gordon equation can also be neglected since 62 1 implies that 

4>/{Hcj)) <C 1. However, in the two-field case, the properties <C 1 and 

(pi/iyHpi) <C 1 cannot be deduced from 62 <C 1. As a consequence, neglecting 
the acceleration in the Klein-Gordon equations for the heavy and light fields is in 
fact an additional assumption that we will make in the following. Then, using that 
H = —K /2 (which, by the way, shows that the Hubble parameter always 

decreases) and the slow-roll Klein-Gordon equation to relate the first time derivative 
of the fields to the derivative of the potential, one obtains 

,2 + 


61 ~ 2Mi 


(R‘4i + <t>j) 


2 • 


(1.189) 


If the heavy field dominates, ii^h 3> <pe or, equivalently, rrih^h m£4>e^ then 
61 ~ 2 Mpj/(/)^. As expected, this expression is similar to that one would obtain in 
single Large Eield Inflation (LEI). And if the light field dominates, i.e. if pi S> R'^ph, 
then 61 ~ 2Mpj/())£. Therefore, we will assume the following initial conditions which 
guarantee that slow-roll is valid 


ph V2Mpi, Pi 3> \/2Mpi, Rph Pi- (1.190) 


The second condition is a priori less obvious so let us discuss it a little bit more. The 
domination of the heavy field comes to an end when Rph = pi- At this transition, 
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the first Hubble flow parameter is given by eit = 2Mp^{l + So inflation does 

not stop provided (fit ^ \/2-\/l + B?Mpi ~ \/2i2Mpi. Since the light field is almost 
constant during the phase dominated by the heavy field (see below), this justifies our 
initial condition. But it is also possible to consider a situation where inflation stops at 
the transition, namely < y/2RMpi [105-107]. After the transition, ei ~ 2MpJ(f?^ 
and if one wants inflation to start again, one needs the condition (1.190) for the light 
field. 

As long as the slow-roll approximation is valid, the equations of motion can be 
integrated and the solution for the field vacuum expectation values reads [105-107] 


/ 45 

(fh = \ — sin [^(s)] , 
V K 

while the Hubble parameter is given by 

2s 


= \ — cos [0(s)] 
V K 


H‘^{s) = —ml [l + (i?^ — l) sin^ 9\ . 

O 


(1.191) 


(1.192) 


This is a parametric representation of the solution in terms of the variable s defined 
by s = — ln(a/aend )5 with 


s 


^°(cos0)^”^h/(™h-”^l)' 


(1.193) 


The initial phase, dominated by the heavy field, corresponds to 0 —>■ 7r/2, s —>■ oo, 
ei —?■ 0 and (fti/cfi = tan0 —?• oo. As already mentioned, this happens when 
> ini4’£ or 6 < 6t where the “transition angle” 9t is given by tan^t = R~^. 
The fact that R > 1 implies that 9t < ir/4:. If 72 ;i> 1 then 9t — R~^ is a small angle. 
In that case (i.e. 72 3> 1), in the regime where 0 S> 0t (he. 6 is not a small angle), 
we have the following behavior for s: s ~ sg cos“^ 9. This implies that the heavy and 
light fields are given by 


(fh — a/—^ tan6< » Mpi, 

Kj 


4so 


K 


(1.194) 


In this regime, the heavy field is super-Planckian and the model is effectively 
equivalent to large held inhation (LFI). This is conhrmed by writing the Friedmann 
equation (1.192) using Eq. (1.194) 


77^ 


2 d 2 ■ 22so 2. 2n ^ 2/2 

—m^R sm 9 = -—m^tan 9 = — 
o o D 




(1.195) 
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which is exactly the Friedmann equation for LFI. On the other hand, as announced 
above, the light field is frozen and its back-reaction is negligible. This provides an 
interpretation for the parameter sq: it is nothing but the vacuum expectation value 
of the frozen light field. Let us also notice that the condition (/>£ 3> \/2Mpi translates 
into So 1/2. Moreover the condition for avoiding an interruption of inflation reads 
So > i?V2- 

Then, the next question is to calculate the behavior of the two scalar fields after 
the transition. The light field now drives the expansion of space-time. The situation 
is a little subtle because one can still have H <C but 3Ff(/h / —w-h'/h- In other 
words, the background still inflates but the heavy field, that has become a test field, 
is not necessarily in slow-roll. In that case, this, in principle, invalidates Eqs. (1.191), 
(1.192) and (1.193) since they all assume H ^ and the two fields in slow-roll: in 
other words, having the kinetic terms negligible in the Friedmann equation and only 
one field in slow-roll is not sufficient to derive Eqs. (1.191), (1.192) and (1.193). In 
that case, we need to return to the exact Klein-Gordon equation for the heavy field. 
If we write /)h = then it takes the form 


A 




/h = 0. 


(1.196) 


Since the background is still in slow-roll, one can neglect the term H in the above 
equation. Then, we see that the behavior of the field depends on the ratio Since 

H is decreasing, the term proportional to the mass necessarily becomes dominant at 
some time and then the field oscillates, namely 


(/h — a cos (mht). 


(1.197) 


The frequency of the oscillations is given by the mass of the field. The amplitude of 
the oscillations decreases as oc and, therefore, the heavy field becomes negligible 

very rapidly. During the oscillations of the heavy field, inflation continues driven by 
the light field. It comes to an end when the vacuum expectation of the light field 
becomes sub-Planckian. 

Having described the behavior of the background, we can now turn to the 
perturbations. They are described by the Bardeen potential already introduced before, 
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and the two perturbed scalar fields 6(j)h and 5(j)£. The corresponding equations of 
motion read 




+ 

s'(j)£ T 3H5(I)£ T 



Hi 


- {4>H4)h + <j)i64>i^ , 
4<f£^ — 2m'j^(j)£, 


A(j)£^ — 2rnj^(j)£. 


(1.198) 

(1.199) 

( 1 . 200 ) 


Unfortunately, this system of equations cannot be solved analytically. However, on 
large scales, namely for wavelengths larger than the Hubble radius, the expression of 
the growing mode of the Bardeen potential and of the two perturbed scalar fields can 
be established. They read [105-107] 


= -Cn-H 






d /4Hh + d,uA 

dt V Vt + V£ J ’ 

/ dhUh + d£V£ 

( U + r, 

( Vh + r, 


( 1 . 201 ) 

( 1 . 202 ) 

(1.203) 


where Ci(k), dh(k) and d£(k) are integration constants. At this point, the following 
remark is in order. We have seen that, in the theory of cosmic inflation, the source 
of the perturbations are the quantum vacuum fluctuations. This of course remains 
true in a model where we have several scalar fields. This means that the quantities 
$, 5H and 5(f)£ are in fact quantum operators. A convenient way to describe this 
situation without introducing all the machinery of quantum field theory is simply 
to write that the amplitude of the perturbed fields at Hubble radius crossing are 
given by 5H = H/y/2kf e\^{k) and 5(f)£ = H/y/2k^e£{k), where eh(fc) and e£{k) are 
two independent Gaussian stochastic processes satisfying (eh(fc)) = {e£{k)) = 0 and 
(eh(fc)eh(fc')) = 6^^\k-k'), {ee{k)e£{k')) = 6^^\k - k'), {eh{k)e£{k')) = 0. This 
parametrization raises in fact non trivial questions such as the quantum-to-classical 
transition of quantum cosmological perturbations but, in this review, we will not 
discuss these issues [116-118]. 

Then, let us simplify the expression of the perturbed heavy scalar field, see 
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Eq. (1.202), by using the explicit form of the potential. One arrives at 


4>h 


Cl 

H 


2H {de — dh) 






H 


2H {d£ — dh) 


™'h'^h + '^£^£ 




^ + 2 HC 3 , 


(1.204) 

(1.205) 

(1.206) 


where C 3 = di^ — di and where, in the last equality, we have assumed that the light held 
was dominant (namely the second phase of inhation). Then, one can use the slow-roll 
relation — —w-h'l^h and obtains 

^ ~ ^C'3"ih</>h ^ -^C'3m^())h, (1.207) 


where, in the last equality, we have used the fact that, before the onset of oscillations, 
H ^ nii^. Then, as already mentioned above, the held starts oscillating and the slow- 
roll approximation is no longer valid. As a consequence, the above equations can no 
longer be used. During the oscillations, one has equipartition between the kinetic and 
potential energy. As a consequence, (ph) — This implies that 6 ph — m^2(^h<^</>h 

and, therefore. 


Sph r,^4>h 

- ~ 2 -. 

Ph 4>h 

But, in fact, the perturbed Klein-Gordon equation for large scales modes, if one neglects 
its right hand side, is the same as the background Klein-Gordon equation provided the 
potential is quadratic in the held (which is precisely the case in the present situation). 
As a consequence, 54>h is in fact always proportional to 4>h, the slow-roll approximation 
being satished or not. In other words, and hence dph/Ph are constant. So if 

we assume that the heavy held decays into cold dark matter after its oscillations, one 



has [105-107] 


<^cdm 


f^/Oh 

_ 

=2^ 

Ph 

end osci 

start osci 


start osci 


-^C3{k)ml. 


(1.209) 


We conclude that, if one is able to express the constant C^{k), then one can establish 
the expression of the cold dark matter density contrast. But this is in fact an easy 
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task. Indeed, reproducing the same calculation as the one which led to Eq. (1.205), 
one obtains 




( 1 . 210 ) 


Then, using this formula and Eq. (1.205), one can eliminate Ci{k) and find an 
expression for C 3 {k). Straightforward manipulations lead to 


C3ik) = 


1 

'm 


54>i 


<^h 


( 1 . 211 ) 


The next step is to replace the derivatives of the fields by their slow-roll expressions and 
by Heh/y/2k^ (and a similar expression for dcff) as was discussed after Eq. (1.203). 
One arrives at 


Einally, one can estimate the entropy perturbation. According to the definitions 
introduced before, see Eq. (1.174), one has 5cdm = >S'cdm 7 = 3 (C dm — C 7 ) = 
dcdm ~ 4(5.y/3 ~ A dm because dpcdm ~ < 5^7 and, during the radiation dominated era, 
/O 7 Pcdm- As a consequence, one has 


5'cdm,fc — 



H [R^(j)^ ^ei{k) - (^j^^eh(fc)] , 


(1.213) 


where the quantities H, and cfi should be viewed here as scale dependent quantities 
since they are expressed at Hubble radius crossing. In fact, as will be discussed in 
Sec. V, their scale dependence permits the calculation of the isocurvature perturbations 
power spectrum. With these equations, one can now predict the CMB temperature 
anisotropies by using Eq. (1.181). But, in fact, the most important conclusion is of 
course that S'cdm.fc 7 ^ 0. This means that, in a model of inflation with more than 
one held, isocurvature perturbations can be produced. This justihes our claim that, if 
non-adiabatic perturbations are observed in the future, a natural explanation will be 
to consider that several scalar helds play a role during inhation. 

In fact, there is even more. Indeed, during the radiation dominated era, the 
adiabatic perturbations can be written as Cao ~ 34>/2 = Ci{k). The constant Ci{k) 
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can also be evaluated easily using the solutions (1.202) and (1.203). The corresponding 
expression reads 

Crd = “f [</>heh(fc) + (l)eee{k )]. (1.214) 

In particular, one has introduced before the power spectrum of the conserved quantity 
Cfc according to (CfciCfcj) ~ 2Tr‘^/kfV{ki)6^^\ki — ^ 2 )^. In the same manner, one can 
define the power spectrum of the non-adiabatic perturbations by 

{kl - k2 ). (1.215) 

But the most important aspect of the above calculations is that adiabatic and 
isocurvature perturbations turn out to be correlated [104, 108]. This means that the 
correlator is non-vanishing. This correlator can be expressed as 

mkAmM) = (fci - fc2). (1.216) 

This is because the expressions of Ck-, see Eq. (1.214), and 5cdm,fc, see Eq. (1.213), both 
depend on Ch and ei. Erom the above dehnition, let us also notice that one can define 
a correlation spectrum by 

- ( 1 - 217 ) 

Let us stress that, when one constrains the amplitude of isocurvature modes using the 
CMB data, it is of course important to take into account the fact that adiabatic and 
isocurvature perturbations can be correlated. As we will see in the next section, this 
was done in the analysis of the Planck data. 

V. INFLATION AFTER PLANCK 

We have previously studied the predictions of inflation for different cosmological 
observables. In this section, we review what is experimentally known about these 
observables and discuss the corresponding implications for cosmic inflation. 

® Notice that this expression is consistent with the dehnition given above, in the text between 
Eq. (1.73) and Eq. (1.74), namely (CfciCfc 2 ) = 2tt^/ kl'P{ki)S^^\ki + k 2 ) because C-te 2 = Cfe 2 - 
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FIG. 4. Temperature anisotropy multipole moments obtained from the Planck 2013 data versus 
the angular scale £ (notice that, for £ < 49, the scale is logarithmic). The gray points denote the 
value of the multipole Ci for each £ while the blue points represent the value of averaged in 
bands of width A£ ~ 31. The red solid line shows the prediction of the best fit six-parameters 
ACDM model. The error bars correspond to ±lcr uncertainties. The lower panel shows the 
residual signal once the best fit model has been subtracted. Figure taken from Ref. [1]. 


The Planck CMB data have been released for the first time in 2013 [1] and, 
more recently, in 2015, new measurements have been made public [4]. Planck 2013 
has measured the CMB temperature anisotropies and the corresponding multipole 
moments CJ'^ are represented in Fig. 4. Let us remind that these quantities are 
defined as follows. After foregrounds subtraction, the Planck measurements can be 
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FIG. 5. Same as Fig. 4 but with the Planck 2015 data. Notice that the quantity Vi is defined 
by Vi = + 1)Ci/{2t:). This plot should be compared to Fig. 4. Figure taken from Ref. [4]. 


used to construct a map of the CMB temperature anisotropy, namely 


5T 


rp (^) — ^ ^ ■) 


(1.218) 


£m 


where are the spherical harmonics and where the vector e specifies a direction in 
the sky. In practice, 5T/T can be expressed as 


6T 


(e) = j 


dfc 

(27r)^/^ 


F{k) + G{k) 


d 

dr]o 


idp,k-e/a{riis^) 


(1.219) 


where dx = a(??iss)ro + a(T/iss) (ho ~ hiss) (^o being Earth’s radial coordinate and r/o 
denoting the present time) is the angular distance to the surface of last scattering and 
the quantity fc/a(hiss) represents the physical wavenumber of the Fourier mode under 
consideration at the time of recombination. The quantities F{k) and G{k) encode 
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the behavior of cosmological perturbations and are called “form factors” in Ref. [119]. 
Already at this stage, we see that the conhguration where the wavelengths of the 
perturbations become equal to the angular distance of the last scattering surface plays 
a preferred role. Then, the two-point correlation function in real space can be written 
as 

.. - 1-00 

=—'y{2£+l)CrPi{cose), (1.220) 

dvr ^ 
e=o 

where 0 is the angle between the two vectors ei and 62 . This expression defines the 
multipole moments . 

The big novelty of the Planck 2015 data [4, 5] is that they not only lead to a more 
accurate measurements of the , see Fig. 5, but they also provide measurements 
of the Fi-mode CMB polarization. One can then define quantities similar to the CJ'^ 
for the correlation between temperature and Fl-mode polarization fluctuations and for 
the Fl-mode power spectrum. The corresponding multipole moments C™ and Cf^ are 
represented in Figs. 6 and 7. 

Before focusing on the consequences of these data for inflation, let us briefly discuss 
their implications for the standard model of Cosmology. It is important to understand 
that the constraints on the cosmological parameters can depend on the model analyzed 
and on the data used to perform the analysis. In 2013, Planck used the temperature 
anisotropy measurement plus the WMAP polarization measurement on large scales 
{i < 23), the corresponding likelihood function being denoted PlanckTT+WP. In 
2015, at least five different likelihoods have been used: PlanckTT utilizes temperature 
data only and is an hybrid, meaning that the temperature likelihood is not the same for 
low multipoles {i < 30) and high multipoles; PlanckTT-blowP makes use of PlanckTT 
and \ow-i polarization data; PlanckTE-|-lowP corresponds to the TE likelihood at 
£ > 30 plus low-^ polarization data only; PlanckTT,TE,EE +lowP makes use of the 
TT, TE and EE likelihoods at £ > 30 and of the temperature and polarization data at 
small scales. Depending on which likelihood is used, the constraints on cosmological 
parameters can slightly change. 

The theoretical framework used to analyze the data is the flat [i.e. /C = 0 in 
Eq. (1.1)1 ACDM model. In order to specify it, we need to know the energy budget 


5T 6T 
^(ei)^(e2) 
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30 500 1000 1500 2000 


FIG. 6. Multipole moments corresponding to the correlation between temperature and FI-mode 
polarization anisotropies obtained from Planck 2015. The red solid line corresponds to the best 
fit obtained with temperature measurements only. The lower panel shows the residual with 
respect to this best fit. Figure taken from Ref. [4]. 

of the Universe, i.e. the photon energy density the neutrino energy density py, 
the baryons energy density pb, the cold dark matter energy density pc and the dark 
energy density p\ (here assumed to be a cosmological constant). Then, in principle, 
on can calculate the behavior of the scale factor a{t) since we know p = ^ • pi in the 
right hand side of the Friedmann equation (1.1). Of course we also need the Hubble 
rate today, Hq or h = Ho/(100km x s“^ x Mpc“^) which is, therefore, another free 
parameter. We also need the power spectrum of scalar fluctuations assumed to be of 
the power-law form P{k) oc Agk'^s~^ where Hg is the amplitude of the fluctuations and 
ns the spectral index. Here, gravitational waves are supposed to be absent, r = 0. So, 
in this simple framework, the perturbations are characterized by two numbers. In the 
case of inflation, we need three parameters, the amplitude of scalar fluctuations and the 









64 


J. Martin: The Observational Status of Cosmic Inflation after Planck 



30 500 1000 1500 2000 

i 


FIG. 7. Same as in Fig. 6 but for the i?-mode power spectrum obtained from Planck 2015. 
Figure taken from Ref. [4]. 

two first slow-roll parameters. We notice (again!) that the parametrization used here 
is different from what we generically obtain from inflation where the power spectrum 
is not of the power-law form and where r is necessarily non-vanishing (but can be very 
small). Finally, we need a parameter describing reionization and we take the optical 
depth r. The interpretation of this parameter is as follows. After recombination, the 
photons are supposed to propagate freely from the surface of last scattering to us. 
However, at the epoch of the formation of the first stars, estimated to be z^e ~ 10, the 
Universe is ionized again. As a consequence, some of the CMB photons scatter off free 
electrons again. The probability to “avoid” this additional scattering is e“’” where 

^^now 

T = a^ TT-edt, (1.221) 

J tre 

is the optical depth. In this expression, cj.p is the Thomson cross-section and rig is 
the number density of free electrons. When additional scattering occur, the direction 












Cosmic Microwave Background Radiation: II JPBCosmo 


65 


of the photon change randomly and this washes out the CMB anisotropy on small 
angular scales, I > ^re, namely ■ For £ < the CMB anisotropies 

are not changed. The value of i^e clearly depends on z^e- The previous considerations 
imply that, on small scales, the amplitude of the fluctuations becomes and 

there is therefore a partial degeneracy between Ag and r. 

In the Planck papers, one of the free parameters is in fact taken to be where 
the subscript “MC” reminds that this quantity is used in COSMOMC. By definition, it is 
equal to be 6-^^ = 100(rg/dA)|approx- Here, is the sound horizon at last scattering, 
namely 




Cgdry, 


where Cg is the sound speed of the baryon-photon fluid, i.e. 


( 1 . 222 ) 


^2 ^Pb—7 

7 ^Ph T 


1 _ 1 1 
3 Ap^ + 3/9b 3 1 + i? ’ 


(1.223) 


where R = ‘ip]^/[Ap^). The quantity dx is, as already mentioned, the angular distance 
to the last scattering surface and naturally appears in the expression of the multipole 
moments, see Eq. (1.219). Therefore, r^/dp^ is in fact the angular size of the sound 
horizon. 6^^ is defined approximately because its value is calculated at a redshift 
which is given by a fitting formula [ 120 ] 


ziss = 1048 


1 + 0.00124 [1 + gi , 


where the function <71 and g 2 can be expressed as 


gi = 0.0783 1 + 39.5 


2\0.763 


1 -1 


52 = 0.560 


1 + 21.1 (Ob/i^) 


2n1.81 


-1 


(1.224) 


(1.225) 

(1.226) 


with lljn = rib + ric- In practice, instead of including h in the list of free parameters, 
we consider 

We conclude that, a priori, we have a 9 parameters: h or p^, Pu-, Ph-, Pc-, PA, 
Hg, ris and r. However, the photon energy density is not a free parameter because it 
is given by 7r^rg/15 where Tg = 2.7255 ± 0.00006 K is the CMB temperature. In the 
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same way, the neutrino energy density is fixed since pu = A^efr(7/8)(4/ll)^/^/9^ with 
A^eff = 3. Moreover, the fact that the spatial sections are assumed to be flat means 
that, say pA, can be deduced from the knowledge of the other parameters. Therefore, 

the “base” model used in the Planck articles is in fact a six-parameter scenario and it 

is sufficient to ht the CMB data. 

Planck 2013 (i.e. PlanckTT-|-WP using the terminology introduced before) found 
the following results (68% confidence limits) [1] 

= 0.022032 ± 0.00028, = 0.1199 ± 0.0027, (1.227) 

lOOe^c = 1.04131 ± 0.00063, r = 0.089l[];[J}^, (1.228) 

ns = 0.9603 ± 0.0073, In (l0^°%g) = 3.089;[J;[j^^. (1.229) 

On the other hand, Planck 2015 with PlanckTT, TE, EE-|-lowP (as already mentioned, 
using the other likelihoods described before would lead to slightly different numbers) 
gives [4] 

= 0.02225 ± 0.00016, = 0.1198 ± 0.0015, (1.230) 

lOOe^c = 1-04077 ± 0.00032, r = 0.079 ± 0.017, (1.231) 

ns = 0.9645 ± 0.0049, In (l0^°%g) = 3.094 ± 0.0049. (1.232) 

The consistency between Planck 2013 and Planck 2015 is evidently very good. 

More involved data analysis can be carried out by opening the parameter space (for 
instance by considering gravitational waves, a running for the scalar power spectrum, 
a time-dependent dark energy equation of state etc ...) and/or adding more data sets. 
In the following, we will describe the corresponding results for the observables that are 
especially relevant for inflation. 

A. Spatial Curvature 

As discussed in Sec. Ill A, see Eq. (1.23), maybe the most important prediction 
of inflation is that our Universe should be spatially flat (although there are contrived 
inflationary models for which this is not true [121]). Therefore, one can follow the 
strategy described above and relax the assumption that the curvature of spacelike 
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sections is flat. Then, the Planck 2013 data plus the WMAP data on large scale 
polarization imply that [1, 2] 


Oyc = -0.058l°;°^6. (1.233) 

If, in addition, Baryonic Acoustic Oscillations (BAO) data are included, one obtains 
Oyc = -0.004 ± 0.0036. 

The Planck 2015 [4] results have confirmed and tightened this conclusion. Indeed, at 
95% confidence level, PlanckTT,TE,EE+lowP leads to Oyc = —0.040lQg4®. If leasing 
data and BAO are taken into account, one arrives at the impressive following result 

= 0.000 ± 0.005. (1.234) 

Therefore, we live in a spatially flat Universe in agreement with one of the most basic 
prediction of inflation. 

As already mentioned, when one relaxes the assumption that the Universe is 
spatially flat, this introduces a new parameter and, therefore, we are no longer in 
the framework of the six-parameters ACDM base model considered before. As a 
consequence, a priori, the constraints on the other parameters may change. This 
is in particular the case of the spectral index res and its significant deviation from 
the scale-invariant case which is very important for inflation. However, Ref. [5] has 
shown that in the framework where Uyc 7^ 0 and where tensor modes are present, the 
constraint on res becomes 


res = 0.969 ± 0.005, (1.235) 

using PlanckTT,TE,EE-|-lowP. The conclusion that the scale invariant case is ruled 
out seems therefore robust. In fact, Ref. [5] has shown that this conclusion is valid for 
other type of extensions such as different relativistic degrees of freedom (the parameter 
Algff defined above, running, dark energy equation of state etc ...). This is of course 
crucial for inflation. 
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B. Isocurvature Perturbations 


Let us now investigate the Planck constraints on isocurvature perturbations [1, 2, 5]. 
We have discussed before, in Sec. IV D, two types of isocurvature perturbations. Firstly, 
there is the effective mode taking into account cold dark matter and baryons entropy 
fluctuations, see Eq. (1.182) denoted, as already mentioned, GDI in the Planck papers. 
Secondly, there is also the Neutrino Density Isocurvature (NDI) mode and the Neutrino 
Velocity Isocurvature (NVI) mode. Each mode is characterized by its power spectrum 
as in Eq. (1.215) and each cross term can also be described by the correlation spectrum 
as in Eq. (1.217). Therefore, the most general situation can be parametrized by the 
4x4 matrix Vab{k) where a = C, Scdi, •S'ndI) 5'nvi with the convention that 
and similar expressions for the diagonal terms. Of course, this matrix is symmetrical. 

Usually, only a 2 x 2 matrix is analyzed and a power law is assumed for each of the 
power spectra with independent spectral index. But this is not the route followed by 
the Planck team. Instead, they have assumed the following phenomenological form for 
Vabik) 


'Pab{k) = exp 



- In A :2 \ 
-Infeay 


\nVab{ki) + 


/ In /c — In /ci 
\ In A :2 — In k 2 


lnPafe(/c2) , 


(1.236) 


where the two scales ki and k 2 are chosen to be /ci = 2 x 10 ^Mpc ^ and k 2 = O.lMpc ^ 
so that the entire Planck window is spanned. The positive definiteness of the matrix 
requires {Vabf < VaaVbb- 

Then, the following quantities are defined 


where 


(f^min; (-n 


(AT)^, (f^min,Cax) 

(^^)tot (^miri) ^max) 


i^'^)ab (^min, 


t=^max 

E 

^ ^min 


{2i + 1 ) 


(1.237) 


(1.238) 


In this expression, the quantity represents the multipole moments calculated with 
the primordial spectrum taken to be Vab- (^^)tot contributions. 

So, in the standard situation, there is just one contribution and the multipole moments 
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FIG. 8. Posterior distributions for the quantities Oab (^min, ^max) introduced in Eq. (1.237), 
inferred from the Planck 2013 data. No statistically significant deviation from adiabaticity is 
found. Notice that TZ = —( is used in this plot. Figure taken from Ref. [2]. 

are computed with = ^C- '^CC ~ means that the perturbations 

are fully adiabatic. 

In Figs. 8 and 9, we have respectively represented the one-dimensional posterior 
distribntion of ccafe (^min, ^max) and the two-dimensional distribution for the power 
spectra Pab(^i, 2 ) for the three modes, GDI, NDI and NVI obtained from Planck 
2013. The conclusion is clear: there is no statistically significant deviation from pure 
adiabaticity. 

In Ref. [5], the constraints on isocurvature modes implied by the Planck 2015 
data have been derived. This work is particularly interesting since one expects 
the polarization data to have a good constraining power on the amplitnde of the 
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f 1 2 ^ 

FIG. 9. Two-dimensional distributions for the quantities 7^ab(fci,2) = 'Pab inferred from the 
Planck 2013 data for a = 72,, GDI (red), NDI (blue) and NVI (green). Again, amplitude of 
isocurvature spectra and correlation spectra are all consistent with adiabaticity. In these plots, 
notice that Vab has indices a = TZ, T with I = GDI, NDI, NVI. Figure taken from Ref. [2]. 

isocurvature modes. In this analysis, uniform priors for and V((^{k 2 ) are 

assumed in the range [l0“®,10“®]. For the power spectrum of the isocurvature 
power spectra, the same choice is made in the range [0,10“®]. Finally, the adiabatic- 
isocurvature correlation function at k = ki is taken in the range [—10“®,10®]. The 
same quantity, but at = ^ 2 , is fixed through an assumption about the correlation 
spectrum, see Eq. (1.217). Ref. [5] restricts itself to scale independent correlation 
spectrum, 

cosA,b = -^^ (1.239) 

y/raa/bb 

in the range [—1,1]. Writing the above equation at k = ki and k = k2 and requiring 
that the value be the same (since the correlation spectrum is scale-independent) allows 
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FIG. 10. Two-dimensional distributions for the quantities 7^ab(fci,2) = 'P^b ^ inferred from the 
Planck 2015 for different choices of likelihoods indicated by different colors (gray, blue and 
red). This plot should be compared to Fig. 9. The six upper plots correspond to a situation 
where we have a mixture of adiabatic (denoted ADI) and CDI modes, the six middle plots to 
a situation where we have ADI and NDI and the bottom six plots to a case where one has ADI 
and NVI. Figure taken from Ref. [-j]. 
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FIG. 11. Posterior distributions for the quantities Ofah (I'min, f'max) introduced in Eq. (1.237), 
inferred from the Planck 2015 data. This plot should be compared to Fig. 8. Figure taken 
from Ref. [5]. 


us to derive the parameter Vab{k 2 )- 

The constraints on Vabiki) and 'Pab(^ 2 ) obtained from Planck 2015 are represented 
in Fig. 10. The constraints on the quantities Q!ab('^min) ^max) are displayed in Fig. 11. 
The conclusions obtained from Planck 2013 are confirmed and even tightened. No 
isocurvature mode is detected and the primordial fluctuations are fully compatible 
with exact adiabaticity. This has of course very important implications for inflation. 
As explained before, this is a non trivial test for single-field slow-roll models. The 
Planck 2013 data were compatible with this simple class of models and did not require 
to introduce additional fields. The results of Planck 2015 do not modify this claim. 
As we are going to see in the next section, this is also the conclusion reached by the 
Planck measurements of Non-Gaussianities. 
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C. Non-Gaussianties 

Let us now turn to the constraints on primordial Non-Gaussianity, see Ref. [3]. 
Before discussing what was measured by the Planck satellite, it is interesting to review 
how the results are sometimes presented in the literature [122], In order to visualize 
the bispectrum, it is convenient to plot the quantity R 7 ^(fcl,/c 2 , fe 3 )(fci^ 2 ^ 3 )^ in terms 
of the ratios x = k^/ki and y = k 2 /ki with the conditions that fci + ^2 + ^3 = 0 and 
ki > k 2 > k^. This immediately implies that 0 < a: < 1 and 0 < y < 1 and, therefore, 
the visualization can be restricted to this square, see Fig. 12. The fact that k 2 > k^ 
means that y > x and, as a consequence, the red hatched region is in fact forbidden. 
Then, since the three vectors fci, k 2 and form a triangle, every edge length is smaller 
than the sum of the length of the two other edges. This means that y > 1 — x and the 
green hatched region is also forbidden. The conditions k 2 < ki + ks (namely y < 1 + x) 
and ks < ki + k 2 (namely y > x — 1) lead to new constraints but outside the square 
[0,1] X [0,1] and, therefore, are not interesting for us. The previous considerations show 
that it is sufficient to plot B'ji{ki, /c 2 , kz){kik 2 k'if‘ in the white, non hatched, region in 
Fig. 12 in order to have a complete representation of the bispectrum. 

It is common practice to single out particular conhgurations. The squeezed triangle 
corresponds to A:i ~ /c 2 k^ which means x ~ 0 and y ~ 1. The equilateral 
conhguration is given by ~ A :2 ~ ^3 or x ~ y ~ 1. The folded case is defined 
by ki ~ 2k2 — 2^3 or x ~ y ~ 1/2. These three configurations correspond to three 
vertices of the white triangle in Fig. 12. Another configuration is the elongated one 
for which A:i ~ A :2 + fca or x + y ~ 1 and is therefore represented by a line in Fig. 12. 
The same is true for the isosceles triangle /si > A ;2 ~ ^^3 or x ~ y. 

Let us now study how the local bispectrum looks like in this representation. From 
Eq. (1.135), one can write 

6 f 1 

Bn{ki,k2,k3) = -^(27r^)2^-^(l + y^ + x^). (1.240) 

0 x'^y 

where V{k) = Ag{k/k^)^^~^ and, for simplicity, we have taken ns = 1 [strictly 
speaking, one should consider the power spectrum of Eq. (1.75) but, in fact, 
this does not change signihcantly the result of this calculation]. It follows, since 
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FIG. 12. Visualization of the bispectrum shape. From the fact that the three vectors fei, 
k 2 and form a triangle, it is possible to faithfully represent the bispectrum in the white 
triangle. Then, different configurations correspond to vertices or edges of that triangle. 

Bn{h, h, k3){hk2k3)‘^ = kfy'^x'^Bn that 

6 1 

Bniki, k2, k3){kik2k3f = -— (1 + + x^). (1.241) 

lU xy 

In Fig. 13 (top left panel) we have represented this function (without the overall factor 
in the above expression). As it is clear from the plot (and also from the analytical 
expression), the local shape peaks at the squeezed triangle. The local shape has been 
constrained by the Planck 2013 data and one obtains [3] = 2.7 ± 5.8 at 68%CL. 

The Planck 2015 data [6] with temperature only implies f^^ = 2.5 ±5.7 and including 
polarization data, one arrives at 


C = 0-8±5, 


(1.242) 


thus tightening the conclusion that the perturbations are Gaussian. 


Another shape that was studied by the Planck team is the equilateral one. It is 
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FIG. 13. Bispectrum for different shape configuration: local (top left panel), equilateral (top 
right panel), orthogonal (bottom left panel). The bottom right panel represents the slow-roll 
prediction computed for a model where V = w?‘4P‘/2. Notice that the absolute normalization 
in these figures is irrelevant. 


defined by 


B'R,{ki,k2,h) 


18 

y 




2 

s 


111 2 

kfk^ k^k^ kfk^ {kik2kz)‘^ 


1 

fci/cyi 


11111 

/ci/cgfcl k2k‘lk\ k2k‘lk\ k^kfk^ /cafelfcf ’ 


(1.243) 


where, again, one has taken tt-s = 1 for simplicity. One can re-express the bispectrum 
in terms of our variables x and y and then multiply by (^ 1 / 02 ^ 3 )^- This gives 

18 1 

B'R,{ki,k 2 ,kz){kik 2 kzf = —f^{2Ti‘^fAl — + I + + 2xy - y - x - y"^ 

5 ^ xy 

—y‘^x — x^ — x^y). 


(1.244) 
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The corresponding bispectrum has been represented in Fig. 13 (top right panel). The 
coefficient has been constrained by Planck 2013 which finds [3] = —42 ± 75, 

a value compatible with zero. With the Planck 2015 data [6] (temperature only), one 
obtains that = —16 ± 70 and, including polarization, 

/^;( = -4±43. (1.245) 


Finally, the last shape studied by Planck is the orthogonal one for which the 
bispectrum can be expressed as 


BTi{ki,k2, ks) 
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(1.246) 


which leads to 


B-jiiki, k2, k-i){kik2k^)‘^ 


^Cc^°{27r^?Al — {3x^ + 3 + 3/ + 8xy - 3y - 3x - 3y^ 
0 xy 

-3y‘^x-3x‘^ -3x‘^y), (1.247) 


and is plotted in Fig. 13 (bottom left panel). The coefficient has been measured 

by Planck 2013 and the result reads [3]: = —25±39. This conclusion is confirmed 

by the Planck 2015 measurements [6], namely = —34 ± 33 (temperature only). 

If polarization data are included, then one finds 


^ortho ^ _ 2 g ^ 21. 

•I IVT, 


(1.248) 


Once again, the measured value is compatible with Gaussian primordial fluctuations. 

It is also also interesting to represent explicitly the slow-roll result using the same 
visualization tools. This bispectrum was derived in Eq. (1.141). Expressed in terms of 















Cosmic Microwave Background Radiation: II JPBCosmo 


77 


X and y, each term reads 
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(1.251) 

(1.252) 


If we write = //^/(16Mpjei)A:f/^*^ = (27r)^ei^gA:^/^*\ then we see from the above 
equations (1.249), (1.250), (1.251) and (1.252) that the functions only depend on 
X and y. In particular, this definition factors out the term k\. As a consequence, using 
Eq. (1.141), the quantity B^{ki,k 2 -,k^){kik 2 k^)‘^ can be written as [89] 


B^^{ki,k2M){kik2k^? 



^ /W(x,y). 

i=l,2,3,7 


(1.253) 


This bispectrum is represented in Fig. 13 (bottom right panel). We notice that it is 
similar (up to a sign) to the equilateral shape (1.243). 

Finally, Planck 2013 has also measured the four point correlation function for the 
local conhguration. The corresponding constrain on the reads [3] 


< 2800, (1.254) 

at 95% confidence level, that is to say, a result compatible with Gaussianity. A 
recent analysis [123] has confirmed this conclusion. Ref. [123] has indeed found 
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= 0.3 ± 0.9 X 10^ and = —1.2 ± 2.8 x 10^. Finally, Planck 2015 [6] obtained 
5 nl ~ (—9.0 ± 7.7) X 10^ at 68% confidence level. 

We conclude this section on Non-Gaussianity measurements as we concluded the 
section on isocurvature modes: the fact that we do not detect a signal beyond the 
vanilla situation is another non-trivial test for single-field slow-roll inflation with a 
minimal kinetic term. In the remaining part of this review, we therefore focus on this 
class of models and derive the corresponding implications that can inferred from the 
Planck data. 


D. Slow-Roll Inflation 


We have seen in Sec. IV A that the power spectra of scalar and tensor perturbations 
can be expressed in terms of the slow-roll parameters. Since the CMB measurements 
constrain the power spectra, they also constrain the slow-roll parameters [124]. 
In Fig. 14, we show the two dimensional marginalized posterior distributions for 
the parameters ei, e2, 63 and P*, where this last quantity represents the overall 
normalization of the power spectrum"^ obtained from the Planck 2013 data. We see that 
P* and 62 are well constrained while there only exists an upper bound on ei and almost 
no constraints on 63. Explicitly, one has 3.035 < In (10^°P*) < 3.15, log(ei) < —2.01 
and 0.023 < 62 < 0.063 at the two sigma level. Planck 2015 [5] has also analyzed 
this question and found ei < 0.0068 and 62 = 0.0291qqq® using PlanckTT-|-lowP and 
restricting the hierarchy at first order in slow-roll. When high-.^ polarization data are 
included in the analysis, one finds ei < 0.0066 and €2 = 0.030^'^gQQg. 

Let us now discuss the physical information on inflation that can be inferred from 
the above results. Firstly, from Eqs. (1.75) and (1.76), one has at next-to-leading order 
on slow-roll 



(1.256) 


That is to say, we have re-written Eq. (1.75) as 



(1.255) 


which defines the quantity P,. 
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FIG. 14. Two dimensional posterior distributions of the parameters ei, e 2 , £3 and P* obtained 
from the Planck 2013 data. The pivot scale is chosen to be = 0.05Mpc“^ and the 

priors are taken to be as follows: a Jeffreys’ prior (i.e. a flat prior on the logarithm of the 
corresponding quantity) for P* such that ln(10^°P») G [2.7,4.2], a Jeffreys’ prior for ei such 
that log(ei) G [—5, —0.7] (the choice of the upper bound ensures that ei < 0.2 and, therefore, 
that the slow-roll approximation is valid) and flat priors for £2 and £3 such that £2 G [—0.2,0.2] 
and £3 G [—0.2,0.2]. Figure taken from Ref. [16]. 


from which we deduce that, at second order in slow-roll, the Hubble parameter during 
inflation can be expressed as [124] 


Ml 


Svr^eiP* [1 + 2(1 -I- C)ei C*e 2 ] • 


(1.257) 


Since we know the posterior of P* and ei, one can derive the corresponding one for 
H*. The result is represented in Fig. 15, see the red dashed curves. Clearly, the fact 
that we only have an upper bound on ei implies that we also only have an upper 
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— PLANCK+BICEP2/log{ei)-prior 


— PLANCK+BICEP2/ei-prior 

-- PLANCK/log{6i}-prior 


-- PLANCK/-prior 




FIG. 15. Marginalized posterior distribution for the inflationary Hubble parameter at the time 
of pivot crossing with a Jeffreys’ prior (left panel) and a flat panel (right panel) on ci (left 
panel). The dashed red line represents the distribution obtained from the Planck 2013 data 
while the solid black line corresponds to the case where the Planck 2013 data are combined with 
the BICEP2 measurement (here, interpretated as a detection of gravity waves) and illustrates 
how a detection of primordial gravitational waves could allow us to determine the energy scale 
of inflation. Figure taken from Ref. [124]. 


bound on H^. With the Jeffreys’ prior on ei (see the left panel in Fig. 15), one obtains 
In ( 10®—< 1.6, that is to say 

< 1.2 X lO^^GeV. (1.258) 

One obtains a similar number if a flat prior on ei is assumed (see the right panel in 
Fig. 15). Those values can be expressed into gravitating energy scales through 

^ 2.2 X lO^^GeV, (1.259) 

where this value assumes a Jeffreys’ prior on ei (again, a similar result is obtained with 
a flat prior). If primordial gravity waves are detected, then this would fix the value of 
r and, hence, the energy scale of inflation. This is illustrated in Fig. 15 where we have 
also plotted the posterior distribution of obtained when the BIGEP2 results [34] 
are taken into account (assuming, for the sake of illustration, that they correspond to 
a detection of gravity waves). 
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Secondly, let us now study what the constraints on the slow-roll parameters mean 
for the shape of the inflation potential [124]. From Eq. (1.34), we see that this gives 
an upper bound on the first derivative of the inflation potential, namely 



(1.260) 


Using PlanckTT+lowP, the recent Planck 2015 data [5] implies that |U^| < 

0.116 (U/Mpi). On the other hand, the second Hubble flow parameter gives infor¬ 
mation about the second derivative of the inflaton potential. From Eq. (1.35), one sees 
that 



(1.261) 


From this expression, we also obtain bounds on the second derivative of the potential. 
Indeed, one has > -e 2 sup and < 2eisup - e2min/2. Explicitly, one 

has 


0.03 < Mp2 ^ < 0.008. 


(1.262) 


Planck 2015 [5], using PlanckTT-|-lowP, finds the following value M^y^^jV = 
—O.Ol^QQQQ at 95% confidence level. 

Thirdly, although the shape of the power spectrum is entirely characterized by 
Eq. (1.75), it is also interesting to derive constraints on the so-called power-law 
parameters [124]. These parameters are in fact simple combinations of the Hubble 
flow parameters, as exemplified by Eqs. (1.86). We now investigate this question in 
more detail. In Eqs. (1.86), we gave the spectral indices at first order in slow-roll. At 
second order, they read 


Us — 1 — (2ei -|- 62 ) — 2e^ — (3 -|- 2 C)eie 2 — 
tI'y = —2ei — 2eJ — 2(1 -|- C)eie2, 


(1.263) 


while the tensor-to-scalar ratio can be expressed as [see also Eq. (1.87)] 


r = 16ei(l -I- (762). 


(1.264) 
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- PLAIMCK+BICEP2 / log{ej )-prior 

-- PLANCK/log{e])-prior 
- - PLANCK+BICEP2 / e, -prior 
■ PLANCK/6j-prior 








FIG. 16. Marginalized posterior distributions for the derived power law parameters n^, r, tit, 
Qfs, clt and /3 t obtained by importance sampling from the distributions of the second order slow- 
roll parameters. The dashed red line and the dotted green lines are the distributions obtained 
from Planck 2013. We have also represented the results obtained by combining Planck 2013 
with BICEP2 (see the solid black line and the dotted dashed blue line). The most striking 
feature which would follow from a detection of gravity waves (here illustrated by including the 
BICEP2 results taken at face value) is of course that r, and therefore tit using the slow-roll 
consistency relations (1.267), would now be measured. Eigure taken from Ref. [124]. 
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One can also define the runnings for scalar and tensor and, in the slow-roll 
approximation, they are second-order quantities and their expressions read 

as = -2eie2 - € 263 , = -2eie2, (1.265) 

Finally, let us mention that the running of the running for the tensor mode is also 
completely specified by the first three Hubble flow functions and is given by 

/3t = - 2 eie 2 (es + 63 ). (1.266) 


One sees that, in general, one has six independent quantities, namely r, ns, Ut, as, ax 
and /3 t- However, the predictions of slow-roll inflation can be expressed in terms of 
three Hubble flow parameters (at least at this order), ei, €2 and £ 3 . This implies that 
all the parameters describing the tensor sector can, in fact, be expressed in terms of 
those characterizing the scalar sector. Explicitly, these so-called consistency relations 
can be expressed as 


- “o: 


ax — 




8 L8 

r 


- -I- (ns — 1) , 

^ + (ns - 1) (1 - ns - 0 -k 


(1.267) 


:Q^s- 


As before, since we know the posterior distributions of the slow-roll parameters for 
the Planck 2013 data, we can infer those of the power-law parameters. They are 
represented in Fig. 16, see the red dashed and dotted green lines. We see that the 
scalar spectral index ns is very well constrained and is around ns ~ 0.96 (the Planck 
2013 value, with WMAP large-angle polarization, reads ns = 0.9603 ± 0.0073). On 
the other hand, we only have an upper bound on the tensor-to-scalar ratio which is of 
course expected since r oc ei. At two sigmas, one obtains log(r) < —0.88 which gives 


r <0.13. 


(1.268) 


Notice that this result is obtained assuming a Jeffeys’ prior on ei. If, instead, a flat 
prior is chosen, one has log(r) < —0.64, leading to r < 0.23. This is because a flat 
prior has the tendency to favor large values of r compared to what is obtained with a 
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Jeffrey’s prior. In Fig. 16, we have also represented the results obtained by combining 
Planck 2013 and BICEP2 assuming that this last signal is due to primordial gravity 
waves. Of course, in that case r is determined and, as a consequence, the tensor 
spectral index is also fixed. It is now known that the BICEP2 signal can be entirely 
explained by dust contamination [125] but, nevertheless, it is interesting to see what 
would be the implications for inflation of a detection of primordial gravity waves. 

Recently, Planck 2015 [5] has also put constraints on r. As usual, these constraints 
depend on the data sets used and on the assumptions made about the theoretical 
frameworks. Here we just quote two numbers. Using PlanckTT,TE,EE+lowP and 
considering that r is the only extra parameter beyond the base ACDM model, one 
obtains 


^*0.002 < 0.1, (1.269) 

at 95% confidence limit. If instead PlanckTT+lowP+WP (we remind that WP means 
the polarization data on large scales measured by WMAP), this number becomes 
^ 0.002 < 0.09. Here, the subscript “0.002” indicates that the pivot scale is taken to be 
0.002 Mpc"^ 

Very recently, a joint analysis by the BICEP2/Keck Array team and the Planck 
collaboration was released [7]. The results are presented in Eig. 17. In Ref. [34], 
BICEP2 announced the detection of primordial gravity waves at a level corresponding 
to a tensor-to-scalar ratio of r ~ 0.16. The reason for this claim can be seen in Eig. 17. 
In this plot, the red solid curve is the signal due to the weak lensing of E-mode that 
produces E-modes on small angular scales. This contribution is necessarily present in 
the standard model of Cosmology and its amplitude can be inferred unambiguously 
once we know the value of the cosmological parameters. The black dots represent 
the signal measured by BICEP2 and Keck Array. As is well visible, in the range 
i ~ [50,120], there is an excess of power with respect to the red solid line and, 
hence, there must be another source of E-modes. BICEP2 interpreted this excess 
as a contribution coming from primordial gravity waves. However, there is another 
known source of contamination: dust. BICEP2 could not measure accurately the dust 
contribution because it operates at a single frequency only. The BICEP2 team therefore 
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FIG. 17. B-mode CMB polarization multipole moments compared to the leasing signal (solid 
red curve). The black dots represent the values of Cf^ obtained from the BICEP2/Keck array 
map. On recognizes the bump in the range £ ~ [50,120] that deviates from the leasing curve 
and that was interpreted as a detection of primordial gravity waves with r ^ 0.16 in Ref. [34]. 
The blue dots correspond to the same multipole moments but after subtraction of the dust 
contribution, estimated from the cross-spectrum with the Planck 353 GHz channel. Glearly, 
the new data are in good agreement with what is expected from leasing. Figure taken from 
Ref. [7]. 


used theoretical models available at that time to remove the dust signal. On the other 
hand, Planck measures the CMB at different frequencies and, as a consequence, can 
estimate with good precision the dust contribution. It is therefore clear that a joint 
analysis between the two teams is the best way to use at the same time the good 
sensitivity of BICEP2/Keck Array and the good control of the dust signal of the 
Planck team. The result of this analysis are the blue dots in Fig. 17. We see that the 
bump has disappeared which means that the excess of power observed was probably 
entirely due to dust contamination and not to primordial gravity waves. The signal is 
now compatible with leasing. The new analysis suggests that the best value of r is now 
r ~ 0.05 but with very low significance and r ~ 0 cannot be excluded. In other words, 
there is no longer a detection of primordial gravity waves. In addition, one obtains a 
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new upper limit which is now r < 0.12 (at 95% confidence limit) instead of r < 0.11 
from the Planck 2013 data. Notice that we obtained before r < 0.13 from the Planck 
2013 data, see Eq. (1.268), and not r < 0.11, but this is just due to some differences 
between our analysis and the Planck one (essentially, different priors). 

In Ref. [124], it was demonstrated that the sets of inflationary models preferred by 
Planck alone and BICEP2 alone are almost disjoint, indicating a clear tension between 
the two data sets. Using a Bayesian measure of compatibility between BICEP2 and 
Planck, it was indeed shown that, for models favored by Planck 2013 the two data sets 
tended to be incompatible, whereas there was a moderate evidence of compatibility 
for the BICEP2 preferred models. This means that the three assumptions (i) slow- 
roll inflation is the correct description of the early Universe (ii) Planck 2013 data 
accurately measure CMB temperature anisotropies and (hi) BICEP2 measurement is 
due to primordial gravity waves are mutually exclusive. In other words, if one has 
the theoretical prejudice that slow-roll inflation did occur in the early Universe, then 
Ref. [124] already proved that the value r ~ 0.16 was likely to be overestimated. In 
some sense, the fact that dust contamination can explain the BICEP2 signal reinforces 
our trust in inflation! 

Let us now turn to the scalar running 0 : 3 . At 95% confidence level, one finds 

- 0.012 < as < 0.006, (1.270) 


that is to say a value consistent with no running. Finally, one notices that the quantities 
ttx and /3 t are well-constrained. It is easy to understand why on the example of a^- 
One has ax = (^’/S)^ + (ns — l)(r/ 8 ), see Eq. (1.267). In this equation ns — 1 is known, 
one can take ns — 1 — —0.04 which means that ax — (r/ 8 )^ — 0.04(r/8). This parabola 
has a minimum at r /8 ~ 0.02 which corresponds to ax — —9 x 10“^. The maximum is 
for r ~ 0.13 and gives ax — —0.0004. We therefore expect —0.0004 < ax < —9 x 10“® 
and which (roughly speaking) explains why the distribution of ax in Fig. 16 is peaked 
(see the red dashed line). 

A last remark is in order at this point. Very often, as already pointed out, the 
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power spectrum is parametrized as 

/ u \ns-l+as/2ln{k/kt)+- 

Vcik) = A, , (1.271) 

and a similar expression for the tensors. Clearly, this is not exactly what inflation 
predicts since not expanding ns — 1 = — 2 ei — 62 (if one works at first order in slow- 
roll) in the above formula means in fact keeping an infinite number of higher order 
corrections which is clearly inconsistent since Uq is determined at a fixed order. Of 
course, since ns — 1 is small, for all practical purposes, this does not impact a lot the 
final results. 


E. Model Comparison 

Let us now turn to model comparison [15, 16, 124, 126]. We would like to determine 
the models of inflation that perform the best being given the current CMB data. From 
a statistical point of view, this question is subtle. Indeed, suppose we have two models: 
A4i characterized by one parameter On and M .2 characterized by two parameters 621 , 
622 - What does it mean to claim that model AIi is better than model AI 2 (or the 
opposite)? Naively, one could compare the likelihoods of the two models for the values 
of the parameters leading to the best fits. But model M 2 has one extra parameter 
and, therefore, one expects this model to automatically improve the fit. Therefore, in 
some sense, it would be “unfair” to claim that A42 is better than Mi since it is “more 
complicated”. Moreover, suppose that only for, say, 621 G [10“^*^, 10“^®] does M 2 lead 
to a good while, a priori, 9 i 2 could vary, in say [—1,1]. Suppose, in addition, that 
this does not happen for Mi, namely that for On in its natural range of variation, the 
fit is always “reasonable”. How do we take into account this wasted parameter space 
for model A42 in our assessment of the respective performance of the two models? 

In order to answer these questions, one recalls that if ^ 2 (^ 215 ^ 22 ) = 
p{D\ 92 i, 022 -,M 2 ) is the likelihood of model M 2 {D represents the data, here we 
have of course CMB data in mind), then the probability of the parameters 021)^22) 
can be expressed as (the Bayes’ theorem) [127] 

P{G 21 i 622 \D,M 2 ) = . .. ^ -^ 2 (^ 21 , ^22)7r(6>2l|Al2)7r(6*22|-Ad2), 


(1.272) 
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where vr represents the prior distributions and <5 is a normalization factor which 
depends on the data and the model. We would like to calculate the probability 
p{A42\D) of model Af 2 and, therefore, we expect a similar equation to hold, namely 

piM2 \D) = -^p{D\M2)tt{M2 ), (1.273) 

p{D) 

where p{D) is a normalization factor depending on the data only and vr encodes our 
a priori information about model A^ 2 - Clearly, Eqs. (1.272) and (1.273) have the 
same structure since they represent two applications of the Bayes’s theorem. To make 
progress we need to know p{D\M. 2 ). But this quantity is in fact easy to calculate since 
f p(02i, (^ 22 l^, M2)dd2id022 = 1, Eq. (1.272) leads to 

£{D\Ai2) = j T2(021; ^22)7r(02l|A^2)7r(022| Af2)d02ld022 

= J p{D\d 2 li 622 , ■M.2)'^{02l\-M.2)T^{(^22\-M.2)dd2ldd22 

= p{D\M2). (1.274) 

Of course the previous considerations apply in general and the quantity £{D\M.i) is 
called the Bayesian evidence of the model M.i and its definition reads [127] 

p{D\Mi) = £ {D\Mi) = J d9ijC ie^j) vr (%|M0 . (1.275) 

The Bayesian evidence is often normalized to a reference model A4ref and one defines 
= T(T)|A4i)/T(Il|A4REF)- In that case, the posterior probability of the model 
Mi (for non-committal model priors) can be re-expressed as 

p{M^\D) = . (1.276) 

j Bj^ef 

In the following, we will give B^bf since this quantity is in one-to-one correspondence 
with the probability of the model Mi. In particular, one sees that p{Mi) > p{Mj), 
namely model Mi is better than Mj (or more probable), if £(Mi) > £(Mj) or, 
equivalently, B^pp > B^ef- 

In order to see why computing the Bayesian evidence answers the questions asked 
before and can give a fair estimate of the performances of a model, let us consider the 
idealized following situation. Let us assume that the likelihood function of model Mi 
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FIG. 18. Ini3^j,p for the all single field slow-roll models with minimal kinetic terms. The 
reference model is taken to be the one where the priors are directly chosen on the Hubble flow 
parameters. Each model is represented by a bar, the length of which is directly proportional 
to In (the numerical value of In being indicated on the same line). A bar on the 
left means that InH^pp < 0 and a bar on the right that lni?(^pp > 0. The color code refers to 
the Schwarz-Cesaro Escalante classihcation [48]. The vertical dotted black line indicates the 
Jeffreys’ categories, see the text for more explanations. Figure taken from Ref. [16]. 
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has width 59ii and that the prior is flat and has width A^n. Since vr is normalized, we 
have vr(0ii) = l/A^n. We also assume that the likelihood is more informative than 
the prior, namely 56ii < A^n. Then, Eq. (1.275) is approximately given by 


8{D\Ml) ~ A,max 


SOii 

Mu 


A similar calculation for At 2 leads to 


£{D\M.2) — A2,max 


5921 5922 
M21 M22 


(1.277) 


(1.278) 


For simplicity, one can take the reference model to be model Afi and, of course, one 
has Brbf = 1. For one finds 

pi2 _ -^2,max Mu 5921 5922 

A,max 59 u M21 M22 ’ 

On this last equation, we see that deciding whether model Afi is better or worst 
than At 2 does not reduce to the comparison of the likelihood function at the best fit, 
>C 2 ,max/Ai,max but that this ratio is corrected by a factor which describes how much 
parameter space has been wasted. So the best model is not the one which has the 
largest but the one which achieves the best compromise between quality of the fit 
and simplicity of the theoretical description. 

As explained before, here, we focus on single-field slow-roll inflationary models (with 
minimal kinetic term) only. At this stage, the strategy is clear: one must evaluate the 
Bayesian evidence of each of these models in order to rank them according to their 
ability to fit the data. This hrst requires to identify all models of this type and this 
was recently done in Encyclopaedia Inflationaris, see Ref. [15]. In this work, about 
200 models have been identified. A model corresponds to a specific choice of potential 
and of priors for its parameters. Two different models can have the same potential but 
different priors. Each model is denoted by an acronym according to the terminology 
introduced in Ref. [16] and, in the present article, we just make use of this convention. 
A detailed justification of the priors chosen can also be found in that reference. In 
Fig. 18, we have represented the Bayesian evidence of the different models (being 
given the Planck 2013 data) by an horizontal bar the length of which is proportional 
to IuRAf) see also the caption of Fig. 18 and Refs. [15, 16, 124, 128]. In order to 
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translate the numerical value of the evidence into strength of belief, we introduce the 
Jeffrey’s scale [127]. If | InS^gpj < 1, then the model is in the “inconclusive zone”, if 
1 < |lni?ppp| < 2.5, then it is in the “weak evidence zone”, if 2.5 < jlniJ^pp] < 5, 
then it is in the “moderate evidence zone” and, finally, if |lni3^pp| < 5, it is in the 
“strong evidence zone”. If the reference model is taken to be the best model, then, 
by definition all InS^pp are negative. In that case, the best models are those in the 
inconclusive zone and those in the “strong evidence zone” can be considered as ruled 
out. 

In Fig. 18, we see that the best Planck 2013 model is KMIII and that 52 models 
end up being in the inconclusive zone, namely: KMIII, ESIy^; Bles, MHIg, Big, 
ESI, BIss, KKLTIs, KMIIv'>o, BEs, ESE, ESI^, KMII, HI, BEs, BEs, RGB, 
RGIi/16, BIph, AI, BIis, MHI, SEIj, SFI, KKLTEtg, BRtg, KKLTI, SBI, RGI, SER, 
PSNIoA, SFEi, PSNIft 2 , PSNIob, PSNEi, PSKRc, LE>o, SFR, ESIp SSBI2, PSNIftg, 
PSNIepA, SSBI4, TWEo, RGIi, SFRs, MHIi, PSNIepB, TWF^^, SBR^^^, LI, SFRi. 
They represent ~ 26% of the models analyzed. We also find that 21% of the models 
are in the “weak evidence zone”, 17% in the “moderate evidence zone” and 34% in the 
“strong evidence zone”. Planck 2013 is therefore able to rule out about one third of 
the inflationary models. Model comparison with Planck 2015 data cannot yet be done 
since the scientific products are not delivered. However, given the consistency of these 
two data sets, we do not expect very different results. 

We have seen before that the “winner” is KMIII which is a string inspired model 
with the following potential [129] 


L((/.) = M* 


1 — a 


( _^V^^e/3E/Mpi)4/3 


(1.280) 


where a and /? are two free parameters. Its heir apparent is ESI^ the potential of 
which is given by [130-132] 




(1.281) 


The fourteenth on the list is the Starobinsky model [9, 133], namely 


(1.282) 
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Actually, all these models being in the “inconclusive zone”, the difference between 
their Bayesian evidence is not significant. This means that, for instance, one should 
view the Starobinsky model as as good as KMIII. In fact, the main common point 
between all these scenarios is that they all possess a “plateau-like” potential, meaning 
that —)• 0 as (?!)—)> oo. We conclude that Planck 2013 has been able to 

constrain the shape of the inflationary potential, a truly remarkable achievement when 
one remembers that inflation can take place at lO^^GeV, and certainly something 
impossible to do in an accelerator! 

We have seen before that, in order to explain the data, we do not need to consider 
models more complicated than single field slow-roll inflation with a minimal kinetic 
term. This does not mean that more complicated models are ruled out (in the 
frequentist point of view) in the sense that, with a carefully chosen set of parameters, 
they can lead to good fits. However, from the previous considerations, we see that 
those models must have a very “bad” Bayesian evidence. Computing the evidence of 
those more complicated models is certainly a difficult task (for instance, for models 
predicting a non negligible level of Non-Gaussianities, one would need to take into 
account the higher order correlation functions). But, in fact, we do not need to carry 
out such a calculation which, at this stage appears to be useless. Indeed we know in 
advance that they are much “worse” than single-field models because of their huge 
wasted parameter space [134]. It is sufficient to know that they all are in the “strong 
evidence zone” and, clearly, we are not much interested in knowing the ranking in this 
Jeffreys category since the models are ruled out (in the Bayesian sense) anyway. 

We have seen how the Bayesian evidence allows us to rank the various inflationary 
models. However, two models with a different number of parameters can have the same 
evidence if the extra parameters are not constrained by the data. This is certainly not 
a desirable property as the model with less parameters is clearly simpler and, therefore, 
should be favored. In order to break this degeneracy, we now introduce the Bayesian 
complexity [135]. For a model Alj, it is defined by [135] 

Cl = (-21og£ (%)) + 21og£ {e^^) , (1.283) 

where (•) means averaging over the posteriors and represents the maximum 
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FIG. 19. Inflationary models in the space (A^uc, In S^ef) • Each model is represented by a 
circle (the radius of which has no meaning) with its acronym, taken from Ref. [10], written 
inside. The four panels corresponds to successive zooms towards the best region (indicated 
by the dashed rectangles), the one with 0 < N^c < 1, namely where all the parameters are 
constrained by the data and a large value of the evidence, namely the model achieves a good 
fit without wasting parameter space. Figures taken from Ref. [16]. 


likelihood estimate of the model’s parameters. One can easily show, see for instance 
Ref. [17], that the number of unconstrained parameters, given a data set, can be 
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expressed as 


Kc = ^;aram-Q, (1-284) 

where -/Vparam represents the number of free parameters of model Afj. We see that this 
gives us a new criterion to discriminate the various models since a model such that 
0 < < 1 ought to be preferred. The Bayesian complexities (given the Planck 2013 

data) of all the Encyclopaedia Inflationaris models have been computed in Ref. [16]. 
In Fig. 19, we have represented these scenarios in the space (A^uc, In . It can be 
noticed that, among the models in the Planck 2013 inconclusive zone, those with a 
minimal number of unconstrained parameters are: ESI^, ESI^^/^, HI, Bl 2 s, RGR, 
AI, BIis, MHI, RGI, SERi, LI„>o, SER, ESIp RGIp MHIj, SBR^j„ and SERi. The 
number of preferred models is now 17, that is to say ~ 9% of the total number of 
models analyzed here. Of course, as already remarked, these models are all of the 
plateau shape. The distribution of models in the four Jeffreys categories versus the 
number of unconstrained parameters is summarized in Eig. 20. 

F. Reheating 

We now describe the constraints on reheating that can be inferred from the Planck 
2013 data. This question was recently studied in Ref. [20]. In Sec. IIIC, we have 
seen that, as far as GMB data are concerned, reheating can be entirely described by 
the parameter Rrad, see Eq. (1.71) or, equivalently, i?reh, see Eq. (1.70). For each 
inflationary model. Ref. [20] has calculated the posterior distribution of the parameter 
InRreh- In order to estimate how much reheating is constrained, it is convenient 
to introduce the ratio Avrin /APin . In this formula, A7rin_Rj,gp is the standard 
width of the prior while AV\nR^^^^ is the standard width of the posterior distribution. 
Therefore if AvrinR^^p/A'PinRj.^p = 1, the posterior is as wide as the prior and reheating 
is not constrained at all. If, however, Avrin/APin then the posterior 

distribution is more peaked than the prior and there is information gain. Glearly, the 
larger the ratio Avrinij^.gp/A'Pin^^^p, the more peaked the posterior. 

The prior on InR^eh has to be chosen carefully and must be justified by physical 
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FIG. 20. Occurrences of inflationary models in the four Jeffreys categories for different values 
of the unconstrained number of parameters. Figure taken from Ref. [17]. 


considerations. Clearly, the energy density at the end of reheating must be smaller than 
that at the end of inflation and larger than at the BBN time where p^mc = (lOMeV)^. 
Therefore, we require p^uc < Preh < Pend- For the mean equation of state, we take 
—1/3 < uJreh < 1 since, by definition, reheating is a non accelerated phase of expansion. 
As a consequence, one can show that this leads to 


In 


' 1 / 4 ' 
Pnuc 

iWpi 


< In i^reh < In 


' 1 / 4 ' 
Pnuc 

ikfpi 


+ -lnf^ 

3 I Mpi 


(1.285) 


The order of magnitude of Rreh being unknown, we choose a Jeffreys prior in the 
above range. Notice that this differs from what was done in the Planck 2013 paper [2]. 
Indeed, in that work, specific reheating scenarios were considered such as instantaneous 
reheating or “restrictive reheating” where, apparently without a strong justification, 
the reheating energy density is fixed to lO^GeV. Moreover, it seems that a prior on 
the quantity AAI* was chosen which is clearly awkward since it does not necessarily 
guarantee that the two physical conditions on preh and w^eh discussed previously are 
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Displayed Models: 170/193 





I Q ST, Q ST, , Q ST, Q ST,,, 0 ST, Q ST, .,,31 


FIG. 21. The quantity A7ri„ , quantifying how much the reheating is constrained, 

versus the Bayesian evidence for Encyclopaedia Inflationaris models (each model is represented 
by a circle the size of which has no meaning). The inset shows the posterior distribution of 
the reheating parameter for the ten best Planck 2013 models. Figure taken from Ref. [20]. 

valid^. In the Planck 2015 paper [5], it seems that this weird approach has been given 
up. The new method now seems closer to what is done in the present article. Notice, 
however, that, if the prior on the reheating energy density appears reasonable, only 
specific values of Wreh are considered which is, of course, not the most general case. 
Let us also remark that Ref. [5] introduces an equation of state parameter during 
reheating, denoted iCint, called the “effective equation of state” but without defining it 
precisely. In particular, it is difficult to know if it is equal to the parameter introduced 
in Eq. (1.64), which is the correct parameter that ought to be used and was introduced 

® An additional problem comes from the description of Ref. [10] by Ref. [2]. Indeed it is claimed in 
this last paper that, for large field models where ~ f", Ref. [19] considered only scenarios of 
reheating for which Wreh = (n — 2)/(n + 2), a wrong claim as can be checked directly by reading 
Ref. [19]. 
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FIG. 22. Same as in Fig. 21 but assuming the mean equation of state during reheating is known. 
The prior-to-posterior width for the reheating energy density ln(pi.eh/.^pi) is represented 
assuming four values of the mean equation of state wJreh, namely TCreh = —0.3 (top left panel), 
threh = —0.2 (top right panel), FJreh = 0 (bottom left panel) and w^eh = 0.2 (bottom right 
panel). Figures taken from Ref. [20]. 


for the first time in Ref. [19]. 

In Fig. 21, we have represented each Encyclopaedia Inflationaris model in the 
space (A7rin/^^^j^/A'Pin^^^j^,lni?[^gp): good models are on the right and models 
for which reheating is constrained are on the top. The horizontal dashed line 
^'^inRreh/^^ini?reh ~ ^ locutes the models for which reheating is not constrained. 
In order to globally assess the value of the constraints, we can define the following 
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quantity 





(1.286) 




which is the mean value of Avrin/APin weighted by the Bayesian evidence, i.e. 
the mean value in the space of models. This is a fair estimate since disfavored models 
will not contribute a lot to this quantity due to their small evidence. Numerically, the 
Planck 2013 data are such that 


\ APlni?,eh 



(1.287) 


which, therefore, indicates that reheating is indeed constrained. 

It is also interesting to assume that the mean equation of state is known. In that 
case, the parameter Tireh only depends on the energy density at the end of reheating or, 
equivalently, on the reheating temperature. In Fig. 22, we have represented the similar 
quantities as in Fig. 21 for four different values of uJreh, namely Wreh = —0.3, —0.2,0,0.2. 
The most striking feature of this plot is that, for positive values of Wrehj the models 
tend to cluster around the horizontal line Avrin /APin = 1. This indicates 

that, for those values of the mean equation of state, reheating is not constrained. As 
a consequence, the number obtained in Eq. (1.287) comes in fact from a region in 
parameter space where TCreh < 0. This conclusion makes sense since, for uJreh < 0, 
the dispersion of the predictions in the (r, rig) space is much bigger than for positive 
equation of state. More details can be found in Ref. [20], in particular concrete bounds 
on the reheating temperature for different models. 

Concluding, the reheating phase is already constrained by the Planck 2013 data. 
The precise values of the allowed reheating temperatures depend on the model under 
consideration and on the mean equation of state. If Wreh > 0, the constraints are very 
mild. It is also worth noticing that two identical models with two different reheating 
histories can have different Bayesian evidence. This means that, given the accuracy of 
the CMB measurements, reheating now needs to be properly included in data analysis. 







Cosmic Microwave Background Radiation: II JPBCosmo 


99 


VI. CONCLUSION 

In this last section, we briefly summarize what we have learned about inflation in 
the recent years. Inflation is a “violent” phenomenon since it could occur at energies 
as high as the Grand Unified Theory scale, i.e. ~ 10^® GeV. It is thus quite remarkable 
to be able to say something about physics at such a high energy scale. The picture 
that seems to emerge from the recent high accuracy astrophysical measurements is 
that inflation is realized in its simplest version, namely single-field slow-roll with a 
minimal kinetic term. Additional features, such as the presence of several fields or 
non-minimal kinetic term, which may appear as (natural) consequences of embedding 
inflation in high energy physics, do not seem to be relevant. If, indeed, inflation is 
really realized in its vanilla version, an important challenge will be (is) to understand, 
from the high energy point of view, why these extra ingredients are in fact not present. 
Also, important questions such as the physical nature of the inflaton field remains 
unanswered. 

The shape of the potential is also constrained and appears to be of the “plateau 
shape”, a typical example of this class of scenarios being the Starobinsky model. 
Popular models such as monomial potentials are now disfavored. 

Interestingly enough, inflationary reheating is also constrained by the Planck data. 
The constraints are model dependent and correspond to an average reduction of the 
prior-to-posterior of about 40%. 

Given this situation, what should be done to increase our knowledge of inflation? 
It is clear that in order to measure more precisely the shape of the potential, one needs 
to constrain the values of the Hubble flow parameters e^. So far, we only have a good 
measurement of the scalar spectral index which is a specific combination of ei and 62 , 
namely ns = 1—2ei—£ 2 - To measure ei and separately, one needs another observable. 
A more accurate measurement of the scalar power spectrum cannot really do the job 
since it involves an additional parameter, € 3 , see Eq. (1.265). We are therefore left with 
either the tensor-to-scalar ratio r, which is directly proportional to ei, see Eq. (1.87), 
or the bispectrum which depends on ei and 62 in a different combination than the 
spectral index, see Eq. (1.253) and Eqs. (1.249), (1.250), (1.251) and (1.252). 



100 


J. Martin: The Observational Status of Cosmic Inflation after Planck 


Measuring primordial Non-Gaussianities has one great advantage: we already know 
in advance where one should find the signal. If one dares an analogy, it is like searching 
for the Higgs boson. We know that if it is not found in a specific window, the 
consequences would be drastic. However, the shortcomings is that the amplitude of 
the signal, f^^ ~ 0 . 01 , is so small that it is not clear whether it is technologically 
feasible. On the other hand, improving the limits on Non-Gaussianities could be very 
rewarding. Many non-vanilla scenarios predict 1 reaching this limit could 

allow us to rule out single field slow-roll models! 

Measuring the tensor-to-scalar ratio is the other possibility. It requires to measure 
the tensor contribution which can be done through a detection of H-mode GMB 
polarization. At the moment, there are considerable efforts in this direction. The first 
claim of a detection of primordial gravity waves was of course made by the BIGEP2 
team [34]. The signal corresponds to a tensor-to-scalar ratio of r ~ 0.16. However, as 
later shown by the Planck team and already discussed before, the signal can probably 
be entirely explained by dust emission [125]. Other ground based experiments are 
currently operating in Antarctica such as BIGEP3 & Keck (three channels: 100 GHz, 
150 GHz and 200 GHz, sky coverage of 1 — 2% and resolution of 30'), SPTP 0 I/SPT 3 G 
(90 GHz and 150 GHz, 6 %, 1.2'), in Ghile such as Atacama B-mode Search (ABS) 
(145 GHz, 2%, 30'), Atacama Cosmology Telescope (ACTPol)/AdvACT (30 GHz, 40 
GHz, 90 GHz, 150 GHz and 230 GHz, 6 %, 1.4'), POLARBEAR/SIMONS (90 GHz, 
150 GHz and 220 GHz, 6 %, 3 — 5') and in the Canary islands such as QUIJOTE 
(11 — 20 GHz and 30 GHz, 65%, 15' — 55'). Soon (2016) in Chile, the experiment 
Cosmology Large Scale Surveyor (CLASS) (40 GHz, 90 GHz and 150 GHz, 70%) will 
start taking data. There are also balloon borne experiments such as EBEX (150 GHz, 
250 GHz and 410 GHz, 8 %, 10') and SPIDER (90 GHz, 150 GHz and 280 GHz, 8 %, 
30' — 40') which are operating in Antarctica and Primordial Inflation Polarization 
Explorer (PIPER) (200 GHz, 270 GHz, 350 GHz and 600 GHz, 70%, 10' — 20') which 
will be starting in 2016 in Palestine in USA (Texas). The most efficient of these 
experiments will reach a level corresponding to r ~ 0.01 in the following 3 — 5 years. 
If one wants to go further, one needs space missions. Two projects appear to be 
particularly promising: the Lite satellite for the studies of H-modes polarization and 
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Inflation from cosmic background Radiation Detection (LiteBIRD) [136] selected as 
one of the prioritized projects in the master plan 2014 by the Science Council of Japan 
and the Cosmic Origins Explorer (COrE+) [137] which is a proposal for European 
Space Agency (ESA) M4 space mission. LiteBIRD has a polarization sensitivity of 
~ 4.5^K X arcmin, a resolution of 0fwhm = 38.5' and a sky coverage of 70%. COrE+ 
can be “light” with a sensitivity of ~ 2.5/iK x arcmin and a resolution of 0fwhm = 6' or 
“extended” with a sensitivity of ~ 1.5/rK x arcmin and a resolution of 0fwhm = 4' (in 
both cases, the sky coverage is 70%). With these space missions, one should be able 
to gain one order of magnitude on r and reach r ~ 10“^ in the next decade, assuming 
no delensing. With delensing, one might be able to probe even smaller values of r. 

Using the above analogy, measuring r is like searching for super-symmetry. We 
do not know at which level it should show up (we do not know the super-symmetry 
breaking scale) but it could be around the corner and, hence, technologically realistic. 
In fact, a determination of r would immediately lead to the inflaton field excursion. An 
excursion which is just Planckian corresponds to a tensor-to-scalar ratio of r ~ 10“^ 
that is to say precisely the limit reached by future space missions. Therefore, given 
that r = 16ei, if ei > 10“^, then one should be able to measure it in the next decade. 

Of course, a detection of primordial gravity waves would also impact model 
comparison. It was recently shown in Ref. [138] that this could allow us to rule out 
almost three-quarters of the inflationary models compared to one-third for Planck 2013. 

In conclusion, detecting B-mode CMB polarization and, hence, primordial gravity 
waves, is probably the next challenge for primordial Cosmology. An additional step 
would then be to check the consistency relation, r = —Sht, which would constitute the 
final proof that vanilla inflation occurred in the early Universe. However, if r is very 
small, this measurement might be too difficult. In any case, at the time of writing, 
detecting primordial gravity waves appears to be the next frontier for inflation. Only 
time will tell whether this is true or not. 
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